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Problems  assoc  U  with  lie  vc-f'ical  r  isc-  and  passage-  through  a  water  ■ 
air  interface  of  small  buoyancy  *  propelled  bodies  are  considered  In  genmal  the  methods 
of  analysis  are  approximate,  and  sometimes  even  crude,  but  the  results  are  considered 
indicative  of  what  may  occur  CaU  ulations  are  I  muted  to  c  onsideration  of  prolate 
ellipsoids  and  usually  to  those  of  fineness  ratio  equal  to  4  The  terminal  velocity  and 
length  of  rtse  distance  required  to  approach  jr\  closely  are  calculated  in  Chapters  ii 
and  ill-  Comparison  with  experiment  indicates  fan  agreement  It  is  further  shown 
that  above  a  diameter  of  2  in  the  Froude  number  no  longer  increases  appreciably  with 
body  size.  A  rough  analysis  of  the  added-mass  variation  cf  ellipsoidal  bodies  as  they 
pass  through  the  free  surface  is  devdoped.in  Chapter  IV.  This  Later  esio  -aie  is  used 
to  indicate  the  expected  velocity  and  acceleration  changes  i  f  oiting  bodies,  in  Chapter  V 
These  are  shown  to  diverge  somewhat  from  those  found  experimentally  Besides  the* 
eftect  of  added-mass  variatiotvbuoyancv  and  drag  changes  as  the  body  exits  are  also 
shown  to  have  an  important  effect  on  the  velocity-distance  and  acceleration-distance 
profiles  Finally,  m  Chapter  VI  the  effect  of  entrained  ware', as  has  been  experimentally 
observed,is  considered  Calculations  of  the  displacement  volume  of  the  boundary  layer 
on  the  body  are  shown  to  be  in  fair  agreement  with  the  measured  entrained-water 
volumes.  This  permits  some  consideration  of  this  entrainment  effect  at  other  scaler 
of  water  exit  occurrence  -  the  effects  of  Reynolds  number  and  transition  from  lamina' 
to  turbulent  boundary- layer  flow  arc-  shown  ro  be  significant  Consideration  of  this 
estimate  of  the  entrained  water  permits  a  'evision  in  the  calculated  velocity  and 
acceleration  profiles;  in  Chapter  VI  this  is  shown  tc  lead  to  close'  agreement  with 
observation.  The  entrained  water  is  thus  found  to  have  an  important  effect  on  the 
water  exit  behavior  of  the  bodies,  a  fact  which  fas  been  suspected  since  the  first 
observations  of  this  occurrence. 
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A  study  '>f  w.iicr-cxit  hydrobulltstu.  s  j.s  «i  pti.iscuf  fluid-mri  harm,  s  ri-M  .i 
lias  been  underway  .since-  September  ll>3t>  Since  this  field  is  analytically  difficult  atu 
ri-httively  unexplored  the  basic  approach  has  been  experimental  An  earlier  Mi* 

describes  the  experimental  v.ork  In  association  with  the  experiments,  certain  cal. 
la  (ions  and  tentative  analyses  have  been  performed  to  obtain  indications  of  what  can  p> 
expected  Thus  the  terminal  velocity  and  rise  distance  calculation  of  Chapters  11  and 
III  were  needed  to  plan  the  experiments  Other  calculations  were  made  in  an  attempt 
to  explain  certain  observed  occurrences  and  to  indicate  the  relative  importance  of  th< 
several  factors  entering  into  the  water  exit  phenomenon. 

Although  admittedly  inexact  in  some  respects,  the  analyses  presented  in 
this  report  may  have  value  ;  •.  estimating  expected  occurrences  in  other  situations 
Thus  the  ballistic  analysis  of  Chapter  III  may  find  use  in  other  buoyant  -  rising- body 
problems.  Similarly  the  crude  water-exit  analysis  of  Chapter  IV  and  its  use  in 
Chapter  V,  besides  being  very  suggestive  for  the  vertical  exit  problem,  may  prove 
even  more  useful  in  approximate  analyses  of  the  angle  exit  situation  This  possibility 
is  currently  being  looked  into 

I  1  •  Added  Mass  Effects  of  Exiting  Bodies 

When  a  body  moves  through  a  fluid  it  naturally  causes  a  disturbance.  For 
unsteady  motions  ai  least  part  of  the  disturbance  may  be  treated  as  an  effective  mass 
of  the  body.  The  body  appears  to  have  an  "added  mass"  of  fluid  moving  with  it  In  an 
otherwise  infinite  ideal  fluid,  this  added  mass  is  a  known  function  of  body  size  and 
shape,  although  it  may  be  difficult  to  calculate  (2)  for  a  body  of  arbitrary  shape 

Experimentally  and  analytically  this  study  of  water  exit  has  been  concerned 
with  bodies  shaped  as  prolate  ellipsoids  (plane  ellipsoids  rotated  about  their  major 
axes)  since  the  added  masses  are  well  known  (3)  for  such  shapes.  For  analytical  con¬ 
venience  the  added  mass  of  a  body  moving  in  a  fluid  is  expressed  as  a  coefficient  k 
times  the  mass  of  fluid  displaced  by  the  body;  thus,  added  mass  equals  kpV-(p  is  the 
fluid  density,  -V-  the  volume  of  the  body) 
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Numbers  in  parentheses  indicate  references  as  listed  on  page  63. 
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H.e  problem  of  water  exit  involve*  .1  lxui\  in  .1  *ciim  ■  inf  mile  fluid  Pbe 
1  •  ee  su-faee  limits  the  liquid  medium  A  very  simple  indicatim*  f  whuf  might  be  ex- 
pee  fed  us  a  buoyant  Ixuly  jvisses  through  a  water-air  interface  is  afforded  by  an  elemen¬ 
tary  momentum  analysis  including  the  effect  of  added  mass.  As  indicated  in  big  1 
when  the  body  is  rising  with  terminal  velocity 
in  water  at  some  distance  ( jy j  >  £)  below  the 
surface,  the  forces  acting  on  it  must  be  in 
equilibrium.  1  e.  .  the  buoyancy  is  exactly 
balanceo  by  the  weight  and  drag.  Once  it 
exits,  the  buoyancy  is  lost  and  the  body 
quickly  starts  to  decelerate  due  to  the  action 
of  its  weight  and  the  air  drag  forces  (which 
are  small,  however).  But,  for  short  while, 
it  would  appear  that  the  body  should  have  a 
larger  velocity  in  air  than  it  has  in  water. 

This  is  shown  by  equating  the  momenta  of 
'he  body  and  added  fluid  before  exit  to 

those  after  exit-  D  L 

w  T 

Using  die  subscript  w  to  Fig.  1  Sketch  of  Body  Before  and  After 

Ex  1 1 

’ndtCdte  the  situation  with  rhe  body  in  water 

and  a  for  rhe  body  in  air,  the  momenta  of  the  body  plus  added  mass  M'  are  equated 
in  the  two  situations 

fM  +  M  ’  )  V  =  (M  +  M '  )  V 
w  w  a  a 

It  is  assumed  that  the  momentum  which  may  be  left  in  the  water  in  the  form  of  waves  01 
other  disturbances,  after  the  body  has  exited,  is  relatively  small  and  may  he  neglected 
Using  the  subscript  b  to  designate  the  body,  the  foregoing  momentum 
expression  may  be  written 

fp.  +kp  )-¥-  V  -  (p,  +  k  p  J  V-  V 
Kb  w  u  b  *  a  a 
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where  p  indicates  the  fluid  density.  The  density  of  air  is  negligible  compared  with 
that  of  water  Jf  the  ratio  p  / is  designated  as  the  specific  gravity  S  of  the 
Imdy,  'hen 
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I  fit  .i  |>|vi  r  »-ii  f  l\  ■  '"iilitfii  i  ik  :  ■  i  •  i  i  <1  \  i  I  *  ■  f pi »  iIk  t 1  -d  |»>.  thi  s  .1  n.i )  •,  s  i  s  ;  i.i  s  fx-i  -n  i  i  i  1 1  <! 

i  i  tin  vi  I  i  M\  j  1111,1"  V  i  i  nul  l  iii  in  M  Mir  <  pt.‘i  ’"(I  m.igmi  n  1  •  ■  I  Uu  >  i.ii.  i 

iisu.i  i  I*><l  i-s  i|  (inrni''S  4  ,•  in  I  h  I!  c-i  i  .idilt-d  m.i's  coefficient-  \  (for  longitudiu.i 

lie  M 1 1  ill)  .1 1  *  1 1  I N  J  .mil  1 1  1 1  IN  l.nw  sprulk  gr  I  vifv  f>odi es  a  re  needed  f o  obfa  ; n  a pjJI  ec  .  - 

uble  velm.  1 t\  pimps  I  hus  ('»r  S  >f  H  2.  Ihe  el'  k  i  Mrs  will  nicicasc  try  41  and  14  per- 
>.  enl  i  i  spec  lively 

The  .uiij.il  siiu.iti  m  is  mure  complicated  than  (he  simple  momentum 
nwlvsis  \c  mid  .mply  Ihe  details  of  the  occurrences  as  the  body  passes  through  the 
interface  are  of  prime  interest,  but  it  is  exactly  these  which  are  glossed  over  by  the 
m  tmentum  approach  Information  o  needed  on  the  variation  in  added  inajs  iii  this 
situation  in  order  to  permit  one  to  evaluate  the  velocities  and  accelerations  to  which 
the  body  will  be  subletted 

I  2  Added  Masses  of  Bodies  m  Infinite  Fluid 

The  added  masses  of  simple  bodies  in  an  infinite  ideal  fluid  have  been 
calculated  (1.  4).  For  a  cylindet  moving  perpendicular  to  its  axis  k  is  1,  but  for 
ellipsoids  moving  paiallcl  to  thin  major  axes  k  vanes  front  0.  S  for  the  sphere  to 
)  HISS  for  a  body  of  fineness  ratio  of  111  File  fineness  ratios  f  d  -  f  of  the  bodies 
represent  the  ratio  of  their  lengths  P  to  maximum  diameters  d  The  variation  of 
k  with  f  d  has  been  tabulated  bv  Lamb(NlH);  this  together  with  the  welted  surface 
area  are  shown  as  a  function  of  P  d  in  the  Appendix  (Fig.  Al) 

Although  the  added  mass  can  be  calculated  in  many  cases  or  at  least, 
estimated  (2)  there  has  been  little  experimental  check  on  these  results  File  added 
mass  is  not  measure-able  --  at  most  it  can  be  inferred  from  experiment  About  the 
mly  determinations  have  been  for  vibrating  bodies  Even  the  simple  case  of  a  body 
accelerating  from  rest  has  had  little  experimental  verification  (S)  Some  evidence  along 
this  line  was  obtained  in  the  experimental  phases  of  this  study  and  was  submitted  in  an 
earlier  rejxirt  (I).  where  tins  pniiil  was  linther  discussed. 

1  ] .  Added  Mass  Variation  Near  Surface 

If  the  boundary  to  the  semi  -  infinite  fluid  were  rigid,  its  effect  would  be 
relatively  easily  analyzed  Thus  hydrodynamics  (S,  4)  indicates  that  tor  a  sphere 
moving  pcrpcndicului  to  a  [dune  rigid  surface 

1  ll’> 

—  -j-  ) 
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k  ( )  S  (  I  * 


(2) 
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whe*e  c!  diameter  of  sphere  uiul  y  the  distance  _»f  its  i  ee.tt «.id  ti  m  flu  surface 
(taken  positive  here)  Since  k  increases  as  y  dee » cases  the  effect  is  that  t  A  a 
repulsive  force.  For  motion  parallel  to  the  surface,  'he  added  mass  decreases  and  it 
appears  as  though  theTe  were  an  attractive  fore”.  The  above  tela'ton  for  k  is  no 
longer  correct  even  before  the  sphere  touches  the  surface  iy  -  d/  2,  hut  it  is  still  of 
v  a  ue  in  indicafr.g  the  order  of  the  effects.  Hie  indu  ite''  k  is  then  me r  eased  fry  32 
percent  For  a  light  body  fS  r  0  2)  this  would  result  in  an  11  percent  rtdut  Mi  n  in 
velocity,  according  to  the  velocity  relation  of  Ecj  i. 

A  free  surface,  cr  interface  between  fluids  of  different  density.  Oc.es  not 

usually  behave  as  a  rigid  wall  since  it  is  free  to  deflect  unde-  the  -estraintng  influence 

of  gravity.  It  is  not  surprising  therefore  that  the  mode  of  behavio-  can  ft*  she  wn  to 

depend  upor.  rhe  Froude  number,  the  ratio  of  ine-rial  to  gravity  forces  '.).  It  has  been 

stated  that  at  low  Froude  numbers  the  surface  behaves  as  though  it  wre  rigid  ''velocity 

normal  to  surface  zero)  while  at  high  Froude  numbers  the  surface  p'essure  must  he 

constant  and  rhe  surface  is  free  to  deform  (zero  velocity  potential  on  surface/  For  the 

later  case  the  added  mass  decreases  as  a  body  approaches  normal  to  a  free  surface. 

This  is  opposite  to  the  effect  noted  above,  for  fhe  sphere  approaching  a  rigid  wall.  In 

the  first  report  on  this  project  (!,'  rhe  surface  behavio'  effects  were  me  re  thoroughly 

discussed  and  reference  made  to  rhe  recent  work  of  Motora  <1^53),  Eisenher.g  1.1*435) 

and  Breslin  and  Kaplan  ( i ^57)  It  was  noted  that  a  first-crder  approximation  to  fhe 

high-Froudc  number  situation  is  found  from  the  inverse  hydrodynamic-image  analysis. 

This  leads  to  an  expression  for  the  added  mass  coefficient  differing  only  m  sign  from 

the  rigid  (low  Froude  number)  result.  Thus  for  the  sphere 

o  .3 

k  =  0.  5  M  -  —  )  1 2a  ' 

64  y 

When  the  sphere  just  touches  rhe  free  surface  (v  ~  d/2).  k  is  reduc  ed  It/  50  percent 
•o  a  value  of  0.  25-  This  would  result  in  a  velocity  mciease  of  55  po-cent  ft  -  the  light 
body  (S  -  0  2t  considered  previously. 

I  4  Some  Experimental  Observations 

A  brief  indication  of  the  nature  of  the  observed  behavior  of  the  exiting 
bodies  appears  in  order.  This  represents  only  fhe  highlights  of  the  results  summarized 
in  reference  (Jiand  is  included  to  help  give  this  report  proper  perspective  The  calcu¬ 
lations  of  the  following  chapters  will  he  compared  directly  with  results  from  the-  exper*- 
ments,  where  possible;  in  other  cases  rhe  comparison  will  be  qualitative 


I  Ik  experiments  have  been 
i  or.ilui  led  with  buoyant  bodies  of  dlairi 
melers  from  I  to  4  in.  and  fineness 
ratios  I  f  d  ol  4  and  8.  A  photo¬ 
graph  of  a  body  after  exit  is  presented 
in  I  ig.  -•  I  he  motion  of  the  bodies  Is 
recorded  by  high  speed  photographs. 

Position  measurements  from  the  films 
are  analyzed  to  yield  velocity  and 
acceleration  profiles,  b'lg.  3  presents 
j  sample  set  of  profiles.  The  velocity 
jump  is  seen  to  he  considerably  different 
from  the  prediction  of  the  simple  momentum 
analysis.  I  he  peak  velocity  occurs 
below  the  free  surface  and  is  only  about 
14  percent  of  that  predicted  by  iicj.  1. 
i  he  photograph  of  b  ig.  2  indicates 
another  deviation  from  elementary  ex¬ 
pectations  --  a  considerable  amount  of  Fig.  2.  Photograph  of  Exiting  Body, 

water  is  entrained  by  the  body  and  exits  with  it. 

I  he  analysis  of  the  experimental  studies  reported  (1),  included  an  attempt 
to  evaluate  (lie  added  mass  of  the  bodies  as  they  approached  the  free  surface.  The 
apparent  added  mass  of  the  body  when  its  nose  just  touched  the  free  surface  was  found 
to  he  a  function  of  Fronde  number  and  it  was  possible  to  define  the  low  and  high 
b'roude  number  ranges.  Figure  4  presents  the  results  (1)  in  which  it  appears  that  k 
decreases  monotonic  ally  with  increase  in  Fronde  number.  A  Froude  number  value  of 
about  3,  based  on  body  diameter  (IF  -  V/VgcT)  separates  the  two  regimes. 


I.  S.  1  he  Problem 

flic  proceeding  discussions  and  the  experimental  results  presented  more 
completely  elsewhere  (1)  serve  to  indicate  the  problem  with  which  this  project  in  general, 
and  this  report  in  particular,  is  concerned.  It  is  apparent  that  rather  small  effects  are 
occurring  as  the  bodies  approach  and  pass  through  the  free  surface.  The  momentum 
analysis  is  rather  too  crude  while  exact  analyses  based  on  potential  theory  are  difficult 
i/l  i|  pile  it  ion,  due  to  the  necessary  considerations  of  die  r\  surface  behavior. 
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Fig.  3.  Typic. 
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Fig.  4.  E.xperi 


I'  u  r  I  he  m  m  •  i  e ,  i  i  1 1  ->  i  p  |  >  i  <  *.i  e  M  s.iys  nothing  about  tile  entj  .t  iin-d  watei  w ;  h  appears  I  .  Ik 
<i  significant  leatui’eol  the  mo  hit vines  J'liis  repo'.  the.  suiiuna  i  i/es  calculaM  m-»  oi 
intermediate  validity  which  give  same  indication  nf  the  relative  impotuiice  >f  the 
various  (actors,  even  though  t rea ting  i  erta in  ones  rather  crudely 

Certain  problems  are  indicated  as  being  im|xn't.ini  and  a  ill  be  considered 
in  the  following  sections  ol  this  repot  Till's  are  the  approach  o!  a  buoyant  body  to  its 
terminal  velocity,  the  variation  in  added  mass,  buoyancy  and  drag  as  the  body  [lasses 
through  the  interface,  and  the  entrainment  of  water  with  the  exiting  !*>dy- 

1  b.  A  c  k  no  w !  edgem  ents 

The  work  repu  ted  is  hardly  the  result  of  the  solitary  work  of  the  author 
Thus  his  colleague,  M.  L  Clark  gave  valuable  advice  concerning  some  >\  the  integrals 
of  Chapter  HI  and  in  general  concerning  other  phases  of  the  work  The  exper  imental 
results  referred  to  at  various  points  were  obtained  with  Prolessor  Clark  (1).  The  senes 
calculations  of  Chapter  III  were  made  by  A.  L.  van  Made  and  J.  C  F  Chow  as  student 
assistants.  As  |iart  of  his  work  as  Research  Assistant  on  the  project  Yau  Wu  collaboiaied 
with  the  author  m  the  added  mass  analysis  ol  Chapter  IV  as  well  as  making  some  ol  the 
calculations  of  rise  distance  (Chapter  III'  The  displacement  volume  calculations  of 
Chapter  VI  were  initiated  by  j.  Chow  as  part  of  his  work  as  Research  Assistant. 
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(3) 


where  represents  the  wetted  or  surface  area  of  the  body.  On  this  basis  Millikan's 

laminar  boundary -layer  analysis*  leads  to  a  coefftclent^for  a  body  of  ftneness  ratio 
(  f/d)  of  8  16jg..  wi»  by  the  relation 


1. 67 


(4) 


where  the  Reynolds  number  IR^  is  based  on  the  length  of  the  body.  This  represents 
only  the  skin-friction  contribution  to  the  drag.  Millikan  noted  that  separation  did  not 
occur  for  streamlined  bodies  and  assumed  that  this  drag  thus  represented  the  total 
drag.  For  comparison,  the  flat-plate  relation  for  this  coefficient**  is  given  by 


1.328 


(5) 


For  the  turbulent  boundary  layer  the  coefficient  according  to  Millikan’s  l/7th  power 
law  calculation  of  the  8.  16  fineness  body  is 


0.077 

T7T 


IR 


t 


(6) 


*  Millikan  calculated  the  variation  in  drag  for  an  airship  model,  termed  the  NPL  long 
model,  which  had  been  tested  in  many  tunnels  (cf.  Dryden  and  Kuethe).  His  results 
are  presented  in  terms  of  a.  coefficient  based  on  an  area  represented  by  the  volume 
of  the  model  to  the  two- thirds  power.  Dryden  and  Kuethe  (8)  give  the  linear  dimen¬ 
sions  as  well  as  the  length,  diameter,  and  volume.  Difficulty  was  experienced  by 
the  author  in  ascertaining  the  surface  area  of  this  model.  Different  methods  of 
estimating  the  area  did  not  agree.  The  more  accurate  method  lead  to  a  value  of 
about  1.6  for  the  constant  in  Eq.  4.  A  less  accurate  method  lead  to  a  value  of  1. 69. 
Study  of  the  material  presented  by  Goldstein  (10,  p.  517)  suggests  that  1.67  is  about 
the  proper  value.  Considerable  attention  has  been  devoted  to  the  results  of  Millikan’s 
analysis  of  an  f  =  8.  16  body  in  this  section,  since  this  represents  about  the  only  com¬ 
plete  drag  calculation  available  for  a  body  of  abiut  the  proportions  of  interest.  The 
results  lie  between  the  well-known  flat-plate  drag  results  and  the  completely  - 
empirical  results  being  considered. 

**  For  an  ellipsoid  of  fineness  ratio  4,  the  friction  calculations  of  Chapter  V  yield  a 
constant  of  1.44  for  the  frictional  portions  of  the  drag.  This,  however  is  estimated 
to  be  only  45  percent  of  the  drag,  the  rest  being  pressure  drag  from  separation 
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More  recent  expressions  for  the  flat  plate  drag  cocfliciciit  by  Coles  .)0S4)  jiuI  Larul- 

q 

weber  (1^531  are  slightly  below  Eq  "  at  Reynolds  numbers  below  l() 

In  estimating  the  drag  ot  a  fine  body  at  the  Reynolds  numbers  of  the  order 
of  interest  here,  tire  major  question  is  that  of  the  critical  Reynolds  number  for  transi¬ 
tion  from  laminar  to  turbulent  flow  Plat-plate  data  on  this  transition  indicates  that  it 
will  occur  at  Reynolds  numbers  from  if)'’  to  3  x  10^  ot  larger.  Assuming  that  body 
curvature  has  little  effect,  a  reasonable  value  to  take  for  calculation  is  if/’  Hoerner 
in  converting  the  available  data  on  the  drag  of  fine  bodies  to  the  type  coeffic  lent 
also  adjusted  tlie  data  (by  the  turbulence-factor  method)  to  loi  respond  to  tha’  tor  this 
same  critical  Reynolds  number  The  question  of  pressure-drag  contributions  depending 
upon  the  s-ream lining  of  the  body  's  tail  is  importan’  The  ellipsoids  of  interest  do  not 
have  the  streamlined  tails  of  the  bodies  considered  by  Hoerner  and  Millikan 

The  most  probable  variation  in  the  coefficient  with  Reynolds  number 

and  fineness  ratio  based  on  the  above  considerations  is  presented  ir  Fiy  o.  This  figure 
is  based  on  Hoernet's  Fig-  5  l  <  (13),  with,  a  few  minor  flanges  The  appropr  ta  te 
flat  plate  relations.  Eqs  5  and  ,  are  shewn  together  with  a  transition  curve  between 
them  for  rhe  mixed  flow  occurring  with  a  transition  Reynolds  number  o'  This 

curve  has  been  calculated  according  to  Prand'Ts  mixed  How  approacl  Also  shown  on 
the  ;t  is  Eq  4  fur  the  8  16  fineness- ratio  body  of  Millikan  This  relation  has  led  to 
the  increase  in  rhe  coefficients  for  rhe  finer  bodies  over  the  value’s  suggested  by 
Hoerner.  Otherwise,  Fig  4  follows  his  estimated  curves  rather  c'oselv  As  an  inde¬ 
pendent  check  on  these  estimates  one  experimental  point  is  shown,  as  determined  from 
Jones  (14)  measurement  of  the  drag  of  fineness-ratio  4  ellipse. id.  The  agreement  w’thm 
a  few  percent  is  good  and  suggests  that  a  reasonable  allowance  has  been  made  for 
separation  from  the  blunt  tailed  ellipsoids 


II  3  Terminal  Velocities  of  Rare  Models 

Whcr.  the  buoyancy-propelled  bodies  reac  t  te'mmal  vehn  ity  they  die  in 
equilibrium  and  rhe  buoyant  force  R  is  Ivj  lan<  ed  by  tin-  wcigl  t  W  and  ft  c  drag  1).^. 
at  terminal  in  an  unlimited  fluid  This 


II  I  l-.KMINiAI.  VT.IjUU  P  <>L  SMALL  HI'OYAN'I  BODIkS 


Sim  «  I  IS  tk-s  1 1  able  In  |)CI  I'll  m  the  ex  per  imhiiI.i  I  studies  of  water  e>  1 1 
w  i  III  lx  id  ies  uppi  ■  >.h  lung  i  lie  lire  mii  l.u  e  .il  a  i  •  instant  veluc  I  ty  1 1  was  nee  essa  ry  te 
kiiow  what  le  I  111  lli.il  velni  ities  Im  c.\|)a  l  ( )iie  method  ol  ae  h lev ing  lh is  deratum 
.valid  he  through  expel  inieiii  however  estimates  made  ol  the  use  distance  required  to 
reach  nea  r  -  lemiinal  velocities  indicate  that  distances  of  live  to  twenty-some  feet  would 
>e  required  lot  the  bodies  ol  interest  These  distances  an  not  easily  available  and  the 
•si  bodies  were  not  always  stable  over  such  lengths  of  travel  Therefore,  it  has  been 
necessary  to  c  a  h  ula  te  the  pi  ubanle  let  m  ina  I  velocities  from  drag  information. 

'I  he  drag  ol  bodies  icsults  from  skm-f  riel  ion  and  normal-pressure 
dist r  i bi 1 1 ions  on  lh**  surface  These  depend  greatly  upon  die  nature  of  the  fluw,  i.e., 
tile  bomida  ry  -  la  vet  development  along  the  body  surface  Since  the  boundary  layer 
ilxuit  these  bodies  is  apt  to  be  fkiiily  tu  <  indent, a  nd  transition  conditions  are  not  too 
clear  even  in  simple  flow  cases,  unequivocal  res  Its  are  not  jxrssible  either  experi¬ 
mentally  or  analytically  hollowing  a  consideration  ol  the  skin  friction  for  elongated 
txuiies  the  available  dtag-c  oeffn  lent  indications  are  evaluated  to  permit  estimates  lor 
the  ellipsoidal  Ixulies  of  interest.  This  permits  calculation  of  the  expected  velocities, 
as  a  function  of  body  size-,  fineness  and  specific  gravity  The  few  experimental 
results,  which  it  has  been  |x>ssible  to  obtain,  ate  also  presented  for  comparison  wi'h 
the  predictions 

II  I  Skip  Friction  of  b me  Bodies 

It  has  been  slated  lay  many  authorities  that  the  resistance  or  drag  of  a 
streamlined  body  is  almost  entirely  due  to  skin  friction  ( c f  (b),  p.  15  or  (8),  p  14) 

Since  the  bodies  of  basic  concern  in  this  study  are  tather  fine  ( £,  d  >  4)  it  was  thought 
at  first  that  the  res  i  statu  e  and  hence  terminal  velocity  could  he  calculated  by  using  well 
established  flat-plate  skin-frit  turn  relations  and  that  the  only  question  would  be  the 
Reynolds  number  value  foi  transition  from  laminar  to  turbulent  flow  tn  die  bouidary 
layer.  'I  here  was  a  jxjssi  hi  1 1  ty  that  a  small  additional  term  might  be  needed  to  allow 
for  separation  of  the  flow  neat  tin.-  tail  of  the  body  Calculations  of  the  terminal  rise 
velocities  under  the  above  assumptions  were  made,  hut  the  subsequent  more  rational, 
but  less  precise,  calculations  have  indicated  that  the  velocities  obtained  were  of  the 
irdcr  ol  40  pet  cent  too  high  for  laminar  flow  and  about  10  pet  cent  too  high  for  a  mixed 
laminat-  turbulent  b”'..’nd.i  ry  -  lay  er  Mow 


.11 

A  much  more  exact  approach  to  deter m iiu r ion  of  'he  lx>d,  drag  would  be 
through  the  calculation  of  the  boundu  r  y  -  la  vet  g'owtb  along  the  curved  su*  face  of  the 
body.  The  first  such  calculations  werc  by  Millikan  p)  in  19.12,  although  an  earlier 
"equivalent  flat  plate"  (but  with  a  varying  velocity  outside  of  the  layer)  i  akulahon  was 
made  by  Dryden  and  Kuethe  (H)  Millikan  used  if  e  Karman  momentum-integral  method 
for  determining  the  rate  of  growth  of  the  lamina r  Ixrnndary  layer  will',  the  velocity 
distribution  approximated  by  a  parabolic  profile*  Since  Millikan  used  the  l  "’th  power  - 
law  velocity  distribution  for  the  turbulent  portions  of  the  boundary  layer  Moore  (9)  in 
1935  presented  an  analysis  in  terms  cf  von  Karman's  newly  developed  logarithmic  law 
for  the  velocity  profile  This  lead  to  a  significant  increase  in  the  drag  coefficient  at 
me  higher  Reynolds  numbers.  Further  discussion  of  this  type  of  calculation  and  com¬ 
parison  of  the  predictions  with  wind  tunnel  measurements  on  various  aHs1  ip  hulls  is 
presented  in  Goldstein  (10). 

The  methods  of  drag  calculation  discussed  above  aa-  still  rather  approxi¬ 
mate.  Analysis  methods  for  boundary- lay er  flows  have  been  considerably  advanced  in 
recent  years.  The  1954  paper  by  Weigbardt  (11)  presents  a  simple  method  for  analysis 
of  the  laminar  layer  which  appears  more  neatly  correct  and  the  situation  in  regard  to 
the  turbulent  layer  is  well  summarized  by  Granville  M2>.  Suitable  estimation  for  the 
drags  of  the  test  models  for  the  present  study  probably  should  he  made  through  such 
boundary -layer  calculations  (cf  calculations  cf  Chapter  V’,.  However,  those  are 
rather  lengthy  and  ate  still  app'oxinate  due  to  necessary  assumptions  about  the  transi¬ 
tion  Reynolds  number  and  behavior  of  'he  separated  flow  at  the  mi!  of  the  body.  The 
approach  is  therefore  made  via  empirical  consideration  of  drag-coefficient  results. 

II  2.  Drag  Coefficient  Information 

The  drags  and  terminal  veloc  ities  of  the  test  bodies  have  been  estimated 
with  the  aid  of  Hoerner  s  ana  lysis  (i3)of  the  available  drag-coeffic  ien;  data.  Hoerner, 
however,  has  transferred  the  <  es  jits  mio  a  skin-friction  type  of  coefficient.  This 
permits  one  to  compare  di  redly  with  the  flat-plate  sk  ir-frir'ion  predictions  The 
drag  is  therefore  represented  by  the  'elation 


*  More  up  to  date  laminar  boundary- layer  calculations  based  or  mere  refined 
velocity  profiles  are  possible  following  the  work  of  P^etscl  (Pi,  as  utilized  m 
Chapter  V 
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B  -  W  =  V  y  (1  -  S)  =  I)  <H) 

wheie  -V  represents  the  volume  and  S  the-  specific  gravity  of  the  body,  while  y 
is  the  soecific  weight  of  the  fluid  (water) 

For  the  bare,  i  c  unfinned,  test  bodies  the  drag  D.j.  in  iicj  H  is  cal¬ 
culated  with  Eu  3  and  C  evaluated  from  Fie  3.  This  has  been  dore  for  the  one- 
i  w  * 

and  two- inch  diameter  test  bodies  (ellipsoids  of  revolution)  of  fineness  ratio  of  4  and  H 
The  calculation  procedures  are  presented  in  a  later  section  The  terminal  velocities 
thus  predicted  for  the  bare  models  are  indicated  in  Fig  ba  In  some  cases  two  possible 
terminal  velocities  are  shown,  one  for  the  model  with  transition  from  laminar  to  turbu¬ 
lent  flow  in  the  boundary  layer  at  about  10^,  as  indicated  in  Fig  5,  and  the  other  as¬ 
suming  the  laminar  flow  continues  along  the  model  with  no  transition  The  values 

in  the  second  case  were  obtained  by  extending  the  laminar  lines  in  Fig.  5  parallel  to 
the  flat-plate  line.  For  the  blunter  body  (  f/d  =  4)  this  may  be  a  poor  assumption 
The  possibility  of  lower  critical  Reynolds  numbers  was  also  considered. 
Assuming  a  flat-plate  transition  Reynolds  number  of  5  x  lO"5  and  that  at  high  Reynolds 
numbers  (say  3  x  10  )  the  values  of  C w  are  the  same  as  shown  in  Fig  5,  suitable 
transition  curves  were  faired  in  on  a  similar  plot  (cf  the  upper  plot  of  Fig  7)  In 
this  case  terminal  velocities  were  calculated  only  for  the  finned  bodies. 


11  4  Terminal  Velocities  of  Finned  Models 

Since  the  test  models  are  unstable  as  they  rise  through  the  water  (I),  it 
has  been  necessary  to  add  stabilizing  fins  at  the  tail,  'llie  surface  area  of  these  fins 
is  not  negligible,  hence  the  terminal  velocities  just  presented  are  too  high-  An  attempt 
was  therefore  made  to  estimate  the  additional  drag  resulting  from  the  fins 

The  fins  are  low-aspect-ratio  flat  plates  Since  they  are  thin  then  drag 
is  calculated  as  though  they  were  flat  plates  Since  these  are  one  fifth  to  a  tenth  as 
long  as  the  model,  the  boundary  layer  on  them  should  remain  laminar  and  Eq  5  has 
been  used  to  calculate  the  additional  drag  D^.  due  to  the  fins.  The  total  drag  is  now 
D  -  t  Dj.  where  represents  the  ixire  model  drag  as  previously  estimated 

In  terms  of  the  drag  coefficient,  (cf.  Eq.  3)  this  can  he  expressed  as 
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where  (he  first  term  is  the  bare  model  value  given  by  Fig  S  and  the  second  term  is  the 
increment  due  to  the  fin  Here  A ^  is  the  surface  area  of  one  side  of  a  fin*  and  C^. 
is  the  mean  flat  plate  friction  coefficient  at  the  fin- length  Reynolds  number  A  simple 
version  of  this  equation  that  indicates  more  direc  tly  the  relative  contribution  of  the  fins 
is 
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For  the  range  of  bodies  and  Reynolds  numbers  of  interest,  the  fin  contribution  to  the 
drag  is  17  to  31  percent  of  the  bare  model  drag.  Tins  leads  to  a  7  to  12  percent  de¬ 
crease  in  terminal  velocity. 

Corrected  values  of  C  in  accordance  with  Eq.  9a  have  been  calculated 

for  ft/d  equal  to  4  and  8  bodies  similar  to  the  test  models  These  calculations  are 

based  on  the  constants  HAf/A  of  0.  2S  for  the  i/d  -  4  body  and  0.  13  for  the 

i/d  =  8  body  with  fin  lengths  (for  Reynolds  number  calculations)  of  0.  16  and  0.08  of 

the  model  lengths.  The  values  of  C  thus  obtained  are  indicated  in  Fig.  7  for 

transition  Reynolds  number  of  5  x  1CP  and  10  .  It  should  be  noted  that  in  accord  with 

Eq.  9,  the  drag  coefficient  Cw  is  to  be  used  with  the  body  surface  area. 

The  terminal  velocity  of  several  bodies  similar  to  the  test  models  have 

been  calculated  using  C  from  Fig.  7  The  results  ate  indicated  in  Fig.  6b  as  a 

w 

plotting  of  terminal  velocity  versus  specific  gravity  S  for  bodies  of  two  fine¬ 

ness  ratios  and  two  diameters  The  range  of  velot  tties  of  interest  is  seen  to  be  from 
5  to  20  fps  for  the  1-  and  2-  inch  bodies 


II  5-  Method  of  Velocity  Calculation 

Initially  the  terminal  velocities  were  calculated  by  trial  and  error  processes 
for  given  body  shape,  size  and  specific  gravity  The  results  were  then  checked  by  the 
quicker  indirect  method  starting  with  the  velocity  and  solving  for  the  specific  gravity 
The  results  were  also  checked  by  another  direct  method  based  on  fitting  portions  of  the 
Cw  vs.  IR ^  curves  with  power-law  equations**  such  that  D  =  K  '/"  Such  relations 


*  The  factor  8  thus  appears  since  there  are  four  fins  with  two  sides  each. 

**  Such  laws  are  used  in  exterior  ballistics  and  are  named  after  the  Russian 

Mayevski,  who  found  them  appropriate  for  projectiles  with  various  values  of  K 
and  n  in  different  velociiy  zones. 


IH 

.ire  discussed  in  Chapter  III  Tins  represents  the  simplest  calculation  method  once 
these  power-law  relations  have  been  determined  Since  the  calculations  were  completed 
a  quick  trial  and  error  method,  working  directly  from  the  C^  vs  IR  ^  plots,  has 
been  developed.  This  last  method  also  permits  one  to  visualize  the  effect  of  various 
factors  on  terminal  velocity  and  hence  will  be  Indicated  briefly 

For  the  terminal  situation  Hqs.  i  and  8  are  equated,  whence 

Cw  T"  Aw  = 

In  terms  of  the  relations  between  and  -V-  and  the  body  length  and  diameter,  the 

following  relations  may  be  introduces!  for  ellipsoids 

V-  =  \  d2  ?  =  v  f'3 

A  =  a  l  d 
w 


d,  2 

z(r 


a  =  2  5 


(cf  Fig.  A  1) 


The  preceeding  equation  now  may  be  written 

C  =  ~  V2  -  (1  -  S)  g  d  f2 
w  2v  6 

where  f  =  £/ d  =  fineness  ratio.  Solving  for  the  velocity,  this  yields  the  relation 


For  ellipsoids 

K  0.647 

This  factor  is  nearly  a  constant  (+  l/2  per-.ent)  since  a  varies  little  with  f 

Thus  the  terminal  velocity  is 

VT  -  0.647  ^  jLp I  (10a) 

1  w 

and  the  Froude  number  at  terminal  is 

IF  =  Vj/Jg <T  =  0.647  V(J  *  S)/C  '  (10b) 


These  suggest  an  independence  from  body  fineness  ratio,  but  of  course,  C  is  a 
function  of  f  as  well  as  1R 


CaU  u  la  lion  of  terminal  veliuty.  for  given  f  and  S.  now  proceeds 

quickly.  Hie  drag  coefficient  is  initially  estimated  from  Fig.  i  or  7  and  a  V.j  value- 

found  from  Eq.  Ida  fins  permits  calculation  of  IR  and  a  new  estimate  of  C 

f  w 

The  result  converges  rapidly 


1 1 •  ft  Experiment'  Ver ification 

After  several  years  of  testing  and  the  introduction  of  an  air  gun,to  get  the 
bodies  up  to  speed  faster,  a  modicum  of  experimental  information  has  been  obtained  on 
the  terminal  velocities  of  the  1-  to  4-  in  diameter  test  bodies.  The  results  obtained 
with  the  gun  assisting  the  bodies  up  to  speed  are  felt  to  be  more  definitive  since  velo¬ 
cities  above,  as  well  as  below,  terminal  c,.n  be  obtained. 

Some  of  the  apparent  terminal  velocities  found  experimentally  for  the 
several  test  bodies  were  presented  in  an  earlier  report  (1)  Based  on  the  specific 
gravity  and  size  of  the  body  the  calculated  terminal  velocities  were  found  and  compared 
This  comparison,  together  with  a  few  additional  values  for  the  larger  bodies  is  sum¬ 
marized  in  graphical  form  in  Fig.  ft  In  three  cases,  where  the  calculated  terminal 
velocity  was  more  than  five  percent  greater  than  the  observed,  the  calculation  has  also 


Fig.  8.  Experimental  Verification  of  Calculated  Terminal  Velocities 
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been  made  on  the  basis  of  a  t  r ideal  Reynolds  numbei  of  ^  x  *(>''  This  leads  to  a 
bracketing  of  file  observed  results  as  indicated  by  the  vertical  da  shell  lines.  Apparently 
transition  occurred  between  S  x  1U''>  and  10^  In  the  ease  of  the  Highest  rveyiiulus  no.nhe; 
test  with  die  4- in  diameter  body  the  water  in  the  test  tank  was  allowed  to  rest  at 
least  a  half  hour  before  the  body  was  "shot".  This  settling  period  seems  to  have  per¬ 
mitted  a  higher  transitional  Reynolds  numhei  to  be  reached,  close  to  the  10^  value 
assumed  for  the  majority  of  the  calculations  Only  two  cases  of  the  calculated  velo¬ 
cities  being  more  than  five  percent  below  the  observed,  or  vice  versa,  are  shown 
There  is  no  explanation  for  the  large  divergences  in  these  cases.  Due  to  the  difficulty 
of  assessing  exact  achievement  of  terminal  velocity  in  the  experiments,  the  occurrence 
of  most  of  the  data  within  a  tolerance  of  plus  or  minus  five  pen  ent  must  be  regarded  as 
a  good  check. 

Another  way  of  verifying  the  approach  used  in  the  calculation  of  terminal 
velocity  is  by  calculating  from  the  observed  V  results.  From  Krj .  10b  the 


Fig  9.  Drag  Coefficients  Evaluated  from  Terminal  Velounts  C om pa  >  ed  with 
Predictions.  (Finned  Ellipsoids  of  f  -  4) 
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The  experimental  terminal  velm  ity  'lata  on  d*e  t  4  bodies  f  a  s  0  us  been  reduced  to 
Cw  values  These  are  shown  in  Fig  0  in  comparison  with  the  virvcs  of  Fig.  6  ex¬ 
cept  for  the  2- in  body  with  S  of  O  KI,  'he  agreement  is  very  good  Early  transition 
(at  perhaps  IK  x  d)^)  is  apparent  fo**  some  of  the  tests  with  the  4-  and  4-in  bodies 

Results  tor  responding  n  I'arsition  at  a  higher  Reynolds  number  f^iO^i  are  seen  to 
have  been  achieved  hv  w.n'ing  f<r  disturbances  ip  the  water  to  decay  . 

Although  less  dalu  l  as  S  cr  obtained  with  'lH‘  fine'  body  d  h).  agreement 
comparable  to  that  slawn  in  Fig  has  been  obtained. 


Ii  7  Size  Effe<.  t  or  Vein,  i r y  and  Froude  Number 

In  numc-ous  hyd'ob.olts'u.  s  situations  n  is  desirable  to  know  how  the 
terminal  veloi  ity  of  a  buoyant  r</:\  vanes  with  its  size  For  the  present  problem  it 
was  desirable  to  study  wa’ei-  ev  phenomena  over  the  widest  possible  range  in  Froude 
number  and  pn  fact,'  the  nter  r  clatioosl  ips  between  d,  V^,  and  IF  were  considered 
briefly  earlier  ( i The  ef  fee's  of  diamete'  or  V,^  and  F  arc  simply  evaluated 
from  Fig  n  wi'h  the  aid  ■  '  Hqs  d)a  and  'Ob  't  seems  apprupt  ia'e  to  note  just  what 
should  occur  as  d  is  in<,  teased* 

For  a  body  of  given  shaj  e  0  -  9  >  < 1 ,  and  spec  i f it  gravity  Eq  10a  and 

10b  show  that 

V.  ^  y/d7r” 

I  St 

F  ~  i/v/C 
r  v. 


•  The  significance  ef  these  considerations  becomes  apparent  when  ii  is  noted  that  the 
tests  of  the*  4-  and  4- in  diameter  bodies  we„c  initiated  in  the  hope  that  Froude 
numbers  higher  ;han  dr  sc  i  blamed  witJ-  the  2- in  body  could  be  obtained.  Since  in 
going  from  the  t-to  2  in  size  file  Froude  numbe;  increased  in  the  tatio  1.6;  it  was 
hoped  that  further  size  i tic  tease  would  achieve  some  fuT'the*  int  tease  in  IF’  Actually 
there  was  negligible  gam  w  0  d-e  4- in  bodv  and  about  a  10  percept  gam  with  the 
4 - 1 n  As  will  he  noted  these  limited  in«  r eases  are  predK  tah,<* 


Fig.  10.  Variation  in  Velocity  and  Froudc  Nur'ht:  with  l)i,. mcir  . 
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Fhus  Horn  Fig  3  >r  7  for  (he  f  4  lx>dy^  IF.^.  should  be  nearly  constant  in  the  IF 
range  It/’  to  10  Hus  is  the  range  covered  by  the  light  2-in  .  i-in  and  4-in  bodo  - 
in  the  water  exit  tests.  They  all  achieved  IF^-  values  in  the  range  8  to  9  In  this 
range  V  is,  of  .  mrse.  proportional  to  the  square  root  of  the  diameter. 

As  a  more  explicit  indication  of  the  size  effect,  calculations  of  V  ar-.e 

IF'.,  have  been  made  for  the  f  -  4  body  at  two  specific  gravities  over  a  consider. ib. 

1 

>\*nge  in  sizes  The  results  are  presented  in  Fig  10.  It  is  apparent  that  the  only  w  ■■ 
to  vary  the  F roude  number  at  diameters  above  2  inches  is  via  changes  in  body  dc*nsit\ 
or  via  delay  in  tm  isitton  (increase  in  transition  Reynolds  number  IR  ).  For  sizes  up  t 
1  3  in  the  velocity  increases  roughly  linearly  with  diameter;  above  this  size  it  vane- 
roughly  as  the  .8,  3  power  of  the  size  The  use  of  sand-grain  boundary -layer  transit]  >> 
devices  to  initiate  premature  turbulent  flow  and  hence  produce  lower  terminal  veloi  in- 
on  given  bodies  is  discussed  in  another  report  from  this  project  (15). 

For  diameters  above  about  one  foot,  a  different  power  applies  but  the 
calculations  have  not  been  carried  far  enough  into  the  turbulent  boundary -layer  region 
to  define  this  Possible  effects  due  to  variation  in  transition  Reynolds  number  are  mu  - 
to  be  appreciable  for  body  sizes  between  2  and  10  inches. 


2S 


lil  THh  KISh  D1STANCH-  -  A  BALLISTIC  PROBLEM 


li  is  obvious  th.it  the1  experimental  studies  of  the  water-exit  behavior  of 
buoyancy- propelled  bodies  should  be  conducted  after  the  body  has  risen  sufficient  dis¬ 
tance  to  reach  terminal  velocity  (1).  Since  there  seemed  to  be  no  available  information 
on  the  determination  of  the  required  rise  distance  to  reach  terminal,  it  has  been  neces¬ 
sary  to  develop  a  method  of  analysis  to  be  used  in  making  such  predictions.  This 
problem  in  external  ballistics  has  been  solved  on  the  conventional  assumption  that  the 
drag  at  any  velocity  is  independent  of  the  acceleration  of  the  body.  The  effect  of  this 
assumption  and  that  of  negligible  rate  of  change  of  the  added  mass  of  the  body  is  thougf  t 
to  be  small.  It  has  also  been  assumed  that  the  rising  bodies  are  stable  in  their  motion 


III  1.  Basic  Theory 

To  predict  the  distance  required  for  a  buoyancy-propelled  body  to  reach 
terminal  velocity,  the  relation  between  distance  and  velocity  of  the  rising  body  is  needed. 
This  is  obtained  from  a  relation  for  the  acceleration  due  to  the  forces  acting  on  the  body 
in  the  unsteady  motion.  As  sketched  in  Fig  1,  the  rising  body  is  subjected  to  three 
forces,  buoyancy  pulling  upward  and  weight  and  drag  downward*.  Until  terminal  velo¬ 
city  is  reached  the  buoyant  force  is  the  Larger  and  the  body  is  accelerated  upward. 
Applying  Newton's  second  law  of  motion  to  the  body, 


B  -  W  -  D  - 


W  _  W  dV 
g  g 


(U) 


The  first  two  forces  are  simply  expressed  in  terms  of  the  volume  and  specific  gravity 
S,  hence 


y  (1  -  S)-¥  -  D 


S  dV 

yI  dF 


where  y  is  the  specific  weight  of  the  fluid  medium,  in  this  case  water.  It  is  assumed 
that,  the  drag  D  is  uniquely  associated  with  the  velocity  of  the  body  through  the  water 
and  is  independent  of  acceleration.  Therefore,  a  simple  Mayevski-type  resistance  law 
is  chosen, 

D  -  K  V°  (12) 


* 


In  these  calculations  the  body  is  assumed  sufficiently  far  from  free  or  rigid  surfaces 
so  that  the  forces  which  appear  there  are  negligible. 


lb. 


and  at  terminal  velocity  V.p,  the  acceleration  is  zero  whence 
(B  -  W)  =  Dt  =  K  Vnr 


(12a) 


Inserting  these  expressions  into  Eq.  11  one  obtains  the  relation 


K  (v!J. 


dV 

dt 


or 


i  -  (v/vT)n 


S  dV  =  £  dV 
gO  ‘  $)  dt  g  dt 


(13) 


where  p  =  S/(l  -  S).  This  yields  a  simple  relation  between  velocity  and  time. 

Before  proceeding  with  the  determination  of  the  relation  between  velocity 
and  distance  one  additional  factor  should  be  added  --  the  effect  of  added  mass.  Tins 
added  mass  must  also  be  accelerated  and  it  is  assumed  that  the  rate  of  change  of  this 
added  mass  is  a  second  order  effect.  The  right  side  of  Eq.  13  must  be  multiplied  by 
the  factor 

k 

a  -  1  +  (14) 


unere  the  subscript  s  refers  to  the  submerged  or  infinite  fluid  value  of  k,  Eq.  13 
becomes 


,  /  v  Pa  dv 

{VT  ’  g  dt 


(13a) 


An  expression  for  the  distance  y  from  the  start  of  the  rise  is  obtained 
by  noting  that  V  =  dy/dt,  hence  from  Eq  13a 


%  -  dt  =  £-a 
v  g 


d  V 
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If  the  ratio  V/VT  is  expressed  as  the  variable  x,  this  becomes 


Pa  V 


dy  - 
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x  d  x 


1  -  x 


(IS) 


The  integral  expression  for  the  rise  distance  y  to  a  point  where  the  velocity  is  x 
times  the  terminal  velocity  (x  <  1)  is  thus 
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The  rise  distance  ts  seen  to  vary  with  the  square  of  the  terminal  velocity 


111.  2  Solution  of  Rise  Distance  Equation 

Evaluation  of  the  integral  In  Eq.  15a  In  closed  form  is  possible  for  n 
values  of  1.  2,  and  3.  The  constant  of  Integration  Is  evaluated  from  the  condition 
y  =  0  for  x  ~  0,  and  the  following  results  are  obtained: 

for  n  =  1 

Y  =  -x  -  loge  (1  -  x)  (Ida) 

for  n  =  2 

Y  -  -  j  loge  (1  -  x2)  (16b) 


for  n  =  3 


Y  = 


1 
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?l0*< 


,1  +  x  +  x 

( - 7- 

(x  -  ir 


\  -l.  HP  .  *  1  /2x  +  1  x 
)  +  J3  tan  ( - ) 

-/T 


2/r 


(16c) 


Of  course,  these  all  predict  an  Infinite  distance  to  reach  terminal  velocity  (x  '  1) 

It  is  intuitively  obvious  that  this  should  be  so,  since  the  approach  to  terminal  velocity 
must  be  asymptotic  as  (B  ••  W)  approaches 

The  exponent  n  actually  varies  from  1.3  to  2.5  depending  on  the  Rey¬ 
nolds  number  range  and  fineness  ratio  of  the  body.  By  expanding  the  integrand  of  Eq. 
15a  in  series  form 
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and  integrating,  the  solution  becomes 
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x^  +  2 
pn  +  ! 


(17) 


Of  course  as 


x  approaches  unity  the  convergence  of  this  series  becomes  very  slow 


29. 


The  series  has  been  evaluated  numerically  lor  n  *  1.  S  and  2.5*.  The  results  ,ir<-  in¬ 
dicated  in  Fig.  11  together  with  the  solutions  for  n  =  1,  2.  and  3  from  Eqs  in. 

The  application  of  the  solutions  of  Eq.  15a,  as  presented  in  Fig  ll  '<>  die 

prediction  of  the  rise  distance  is  quite  simple  when  n  is  a  constant  Thus  (nr  ,j  given 

2 

n,  this  figure**  gives  one  the  value  of  the  parameter  Y  -  yg (fia  V  requires!  i<. 

2 

reach  a  certain  fraction  x  of  terminal  velocity.  As  g/ap  V^,  is  determ ina ie, 
once  the  body  and  fluid  are  given,  y  can  be  dlrecfly  calculated 

III.  3.  Calculation  of  Rise  Distance  for  n  Variable 

Unfortunately,  as  is  apparent  from  Fig.  7,  the  drag  coefficient  for  a  bods 
does  not  vary  simply  (Cw  =  KJRm)  with  Reynolds  number  and  the  exponent  tr  and 
hence  n  is  therefore  a  variable.  For  a  given  body  and  fluid,  n  in  Eq.  15  will  have 
one  value  at  low  velocities  and  another  at  high  velocities.  The  problem  must  therefore 
be  solved  in  step-wise  fashion,  by  making  the  elementary  ballistic  assumption  that  n 
is  a  constant  over  a  certain  range  of  velocities.  Once  this  range  is  exceeded  anothei 
value  applies.  With  each  range  of  n  application,  a  different  effective  value  of 
applies. 

The  stepwise  solution  to  the  rise  distance  problem  is  best  illustrated 
through  an  example.  For  clarity's  sake  a  two  step  case  will  be  illustrated.  Thus,  if 
only  body  skin  friction  entered,  would  be  given  by  the  Blasius  type  relation  foi  the 

laminar  range  and  n  would  have  the  value  1.5.  Transition  to  a  turbulent  boundary 
layer  might  be  at  a  critical  Reynolds  number  of  10^.  For  a  considerable  range  of  velo¬ 
cities  above  the  critical  the  exponent  is  then  about  2.  5.  The  body  will  rise  in  accord 
with  the  n  =  1.5  curve  with  an  effective  or  fictitious  terminal  velocity  V^,  which 
it  would  reach  according  to  the  laminar  drag  law  if  there  were  no  transition  to  turbulence. 
When  however,  the  velocity  Vc  corresponding  to  the  critical  Reynolds  number  is 
reached,  the  boundary  layer  will  become  partly  turbulent.  From  this  point  on  the 

n  =  2.  5  curve  will  apply  for  the  rise,  but  with  an  effective  origin  y  as  though  the 

£ 

body  had  been  rising  from  the  start  with  n  =  2.5.  Of  course  in  this  latter  part  of  the 
rise  the  terminal  velocity  10  used  in  the  rise  calcula  tions  is  that  for  the 

(mixed  flow)  new  drag  law  of  n  "  2.5.  The  general  situation  is  depicted  in  Fig.  12. 


*  Numerical  values  obtained  for  the  solution  of  Eqs.  16a,  b,  c,  and  17  are  presented  in 
the  Appendix  (Section  Vill.5). 

**  For  ease  of  calculation  lines  of  constant  x  were  plotted  as  Y  vs.  n.  These 
plots  are  appended  as  Figs.  A2  and  A3  in  Section  VIII .  4. 
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Fig-  12.  Two  Step  Calculations  of  Velocity-Distance  Relation  for  Buoyant  Body 


For  simplicity  of  calculation  the  rise  distance  may  thus  be  written  for  the 
final  distance,  with  which  one  is  really  interested,  as 

y  =  y(x ,  n)  =  y(x.  2.  5)  -  y(x  ,  2. 5)  +  y(x  ,1.5)  (18) 

c2. 5  10 

The  parenthetical  terms  are  symbolic  meaning  that  y  is  to  be  evaluated  from  Fig.  1 1 
(or  Figs.  A2  and  A3)  for  the  particular  value  of  x  =  V/V^  and  n.  In  the  case  of 
the  critical  point,  x^  is  different  for  the  two  exponents,  even  though  the  critical 
velocity  is  the  same  due  to  different  terminal  velocities  applying.  As  may  be  seen 
from  Fig.  12  the  occurrence  of  mixed  flow  shortens  the  rise  distance  when  n^  is 
greater  than  nj.  This  stepwise  solution  is  obviously  extendable  to  cases  when  more 
than  two  drag  laws  apply.  For  each  step  in  which  a  separate  value  of  n  must  be  used, 
the  calculation  of  the  increment  In  y  is  based  on  the  effective  value  of  the  terminal 
velocity  for  the  particular  drag  Law. 

Study  of  the  relation  between  the  drag  coefficient  as  a  function  of  fineness 

ratio  and  Reynolds  number  (cf.  Figs.  5  and  7)  indicates  that  four  or  five  drag  laws  are 

4  7 

needed  to  cover  the  Reynolds  number  range  of  interest  (about  10  to  10  ).  For  the 
various  bodies  of  interest  and  (flat-plate)  transition  Reynolds  numbers  of  5  x  10^  and 
10  ,  these  have  been  evaluated  and  the  several  laws  are  tabulated  in  the  Appendix 
(section  VIII  3).  All  of  the  exponents  are  within  the  range  of  1  to  3,  for  which  the 


n 


unsteady  nmf  on  equation  (Eq.  IS)  haH  been  solved  Since  most  of  the  rise  distance 
occurs  af'er  the  Indy  has  reached  eight- tenths  of  the  terminal  velocity  (c f -  Fig.  12)  a 
three  step  solution  is  adequate  in  most  cases*,  thus 

V  y(x.  n, )  -  yfx  .  n  )  «-  ,(x  ,  n  )  -  y(x  .  n  )  +  y(x  .  n  )  ( 18a) 

1  c21  .)  c22  2  i  Ml  1 

where  n  indicates  the  exponent  for  the  first  drag  law.  etc.,  x  indicates  ihe 

C1  1 

critical  velocity  separating  law  1  and  2  divided  by  the  effective  terminal  velocity 

applicable  to  law  !  at.J  x  indicates  this  same  velocity  divided  by  the  effective 

1  12 

terminal  velocity  for  law  2.  For  x  in  the  first  term  a  value  of  0.  98  or  so  is  chosen 
depending  just  how  close  to  the  terminal  one  desires  to  make  the  calculation 

HI  4  Results  for  Typical  Test  Bodies 

The  rise  distance  as  a  function  of  the  specific  gravity  of  the  body  has  been 
calculated  for  three  typical  test  bodies,  the  I -in.  diameter  finned  model  of  P/d  71  4, 
the  1-in  P/d  -  8  finned  model,  and  the  2- in  jP/d  ~  4  finned  model  (for  F  =  5  x 

r  ^ 

10  and  10  )  The  results  are  shown  in  Fig  13  Except  for  the  1-in  P/d  -  4  body, 
where  only  3  to  5  feet  of  rise  distance  is  required,  it  is  apparent  that  the  necessary  rise 
distance  increases  with  increase  in  S  and  that  5  to  23  feet  are  needed  The  rise 
distance  also  increases  rapidly  with  length  of  the  body  and  the  transitional  Reynolds 
number. 

A  few  scraps  of  information  on  die  rise  distance  are  available  from  the 
experimental  phases  of  this  project  (1).  Thus  over  4  feet  were  required  for  the  1-in 
P/d  -  8  body  (S  =  0.  27).  In  the  case  of  the  2- in.  f/d  -  4  body  about  6  feet  were 
found  to  be  necessary  to  reduce  the  acceleration  to  a  small  value.  According  to  Fig. 

13  some  7  to  9  feet  should  be  required.  In  the  case  of  the  P/d  =  8  body,  8  feet 
should  be  needed  --  twice  the  four  foot  value  noted  to  be  insufficient.  For  a  1-in. 

P/d  -  4  body  it  appeared  that  experimentally  3  to  4  feet  were  required  to  reach  ter¬ 
minal  velocity.  This  body  had  a  specific  gravity  of  about  0.33  whence  Fig.  12  indicates 
a  needed  rise  distance  of  4  feet.  Considering  the  experimental  uncertainties  involved 
in  ascertaining  when  a  body  is  close  to  terminal  velocity,  it  would  appear  that  a  toler¬ 
able  check  on  the  calculations  has  been  obtained. 


*  This  is  a  fortunate  exception  to  the  usual  balltstics  problem  in  which  the  multiplicity 
of  step-wise  calculations  makes  such  an  approach  inconvenient. 
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Rise  Distance 


In  |»)inf.»*  shoA.fi  on  Fig  i  i  wc.c  found  by  the  multiple  step  calculation 
procedure  of  Ivj.  Iha.  P'e  necessity  of  this  complicated  calculation  is  subject  to 
f|uestion  since  in  some  cases  the  extra  terms  had  little  effect.  In  the  case  of  2- in. 
body  the  rise  distances  have  been  calculated  for  every  tenth  change  in  S  over  the 
range  plotted.  Six  of  the  points  were  also  calculated  by  the  more  exact  procedure, 
the  agreement  is  rather  good  in  four  cases  (errors  of  1  to  3  percent).  In  two  cases, 
however,  the  one  step  method  is  not  too  good  since  errors  of  10  and  IS  percent 
occurred  hven  such  errors  may  not  be  unreasonable  when  one  considers  that  the 
definition  of  the  rise  distance  (even  in  terms  of  x  -  0.  9B)  is  physical!)  vague.  In 
many  cases  a  single  step  calculation  probably  suffices. 


IV 


APPROXIMATE  added  MASS  VARIA  non  in  SDKEACh  REGION 


'i  i desired  to  determine  r h «_•  variation  m  tin-  addcd-m  a*»s  cocffn  lent  of 
an  elongated  bod',  is  if  passes  through  the  free  surface  of  a  body  of  water  It  is  obvious 
that  the  virtual  mass  of  the  body  must  deerease  as  it  passes  form  die  denser  to  lighter 
medium  One  c.in  predict  the  deerease  in  added  mass  of  the  lx>dy  and  therefore  the 
velocity  met  ease  U*om  nu  rrentum  analyses)  by  assuming  that  on  exit  all  of  the  momen¬ 
tum  represented  by  ’he  added  mass  in  water  must  be  transferred  to  the  body  As  has 
been  noted,  'his  approach  does  not  appear  to  hm  too  satisfactory 

The  details  of  the  change  in  added  mass  as  the  body  approaches  and  then 
passes  through  the  free  surface  are  also  of  considerable  interest  As  noted  in  Chapter 
I,  theoretical  analyses  of  the  occ  ur  rentes  are  not  simple.  This  present  note  presents 
some  approximate  analyses  of  the  occu1 - cnees  as  the  body  {Kisses  through  the  free 
surface.  Although  the  analyses  are  crude  and  quite  simplified,  it  is  thought  that  this 
approach  is  of  some  value  in  indicating  roughly  what  happens.  The  analysis  is  adapted 
from  a  water  enmy  (  t  impact  analysis  proposed  by  von  Karman  in  Id2d  (16)  for 
estimating  landmg  forces  on  seaplane  hulls. 


p'lg.  14  B<  dy  m  Various  Stages  of  Exit 


i  V  1 .  Analysis  1 

A  body  of  the  shape  considered  is  sketched  in  F  ig  14  as  it  is  exiting  from 
the  water  into  air.  The  body  is  an  ellipsoid  of  revolution,  with  length  f  and  diameter 
d.  The  coordinate  of  (lie  centroid  is  given  by  y,  which  is  positive  when  the  centroid 
is  out  of  the  water  Ihe  radius  of  the  section  of  the  body  hunted  al  the  water  surface  is 
indicated  by  the  symbol  r.  This  has  the  maximum  value  of  d,  2.  The  density  of  the 
water  is  p  and  the  density  of  the  ait  is  neglected. 

It  is  assumed  that  the  water  surface  remains  plane  as  the  body  pjsses 
through  it.  As  y  increases  in  Fig.  14.  file  added  mass  of  the  body  must  decrease 
since  there  is  less  of  it  in  the  water.  Following  von  (Carman's  two-dimensional  approach 
for  motion  in  the  other  direction,  it  is  assumed  that  die  added  mass  is  directly  given  by 
the  radius  r,  and  that  it  does  not  depend  upon  die  manner  m  which  r  varies  with 
y-  The  added  mass  is  therefore  taken  to  be  that  of  a  sphere  of  radius  r,  which  is 
half  water  and  half  air.  Since  the  mass  of  air  is  neglected,  and  the  added  mass  of  a 
sphere  is  one-half  its  displaced  mass  of  fluid,  the  added  mass  of  the  body  is  expressed 
by  the  following  relation 

o  4  3  7  3 

M'  '  4  '  3  7  r  T  p  r  (19) 


The  added  mass  is  usually  expressed  as  a  fraction  of  the  mass  of  a  volume  of  water 
equal  to  the  volume  V-  of  the  body.  The  added  mass  coefficient  is  thus 


M‘  _  2 r  * 

"  ~Ft 


(20) 


The  radius 


r  is  related  to  the  axial  distance 


2 

r 
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y 


by  the  equation 


(2D 


substituting  thus  into  Eq.  20  yields  the  result 


k 


f  <l/4  - 


2  ..2  A .  s 


(22) 


This  gives  ‘he  added  mass  relation  for  the  body  when  ir  is  more  than  half  way  out  of  the 
water. 

Application  of  Eq  22  to  the  case  when  y  is  negative,  i  e.  ,  when  the 
body  is  less  than  halfway  out,  leads  'o  an  increase  in  k  as  the  body  exits.  ThiSjOf 
course,  is  illogical  The  variation  of  k  with  distance  from  the  surface  should  he 
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such  that  at  large  negative  values  of  y  the  added  mass  has  some  constant  value  In 
the  present  approximation  k  should  start  to  decrease  when  the  nose  of  the  body 
touches  the  free  surface  (y  1  -  £/ 2),  should  be  reduced  by  one-half  when  the  body  is 
halfway  out  (y  B  0).  and  should  reach  zero  when  the  tall  is  just  leaving  the  water 
surface  (y  =  + 1/ 2)  The  solution  for  k  should  be  symmetrical  about  toe  value 
obtained  with  y  =  0.  This  value  Is  kQ  =  d/4f.  This  symmetry  could  be  obtained  by 
expanding  Eq.  22  about  this  value  and  letting  the  extra  terms  have  different  signs  for 
the  two  regions,  thus 


where  the  plus  sign  after  the  first  term  goes  with  positive  y  values  and  the  minus 
sign  goes  with  negative  y  values.  Since  the  equation  is  not  to  apply  for  numerical 
values  of  y/|  greater  than  1/2,  the  convergence  is  rapid. 

The  preceding  arbitrary  modification  of  the  solution  to  make  it  symmetrical 
about  the  halfway  point  can  be  shown  to  have  a  rational  basis.  Consider  a  body  less  than 
halfway  out  of  the  water,  as  in  Fig.  14b.  In  terms  of  the  present  approximation  it  has 
an  added  mass  coefficient  of  k'.  When  this  same  body  is  further  out  of  the  water,  as 
in  Fig.  14c,  such  that  its  ntroid  is  as  far  above  the  surface  as  it  was  below  in  the 
previous  case,  its  added  mass  coefficient  is  k".  When  the  body  was  completely  sub¬ 
merged  (Fig.  1 4a)  its  added  mass  coefficient  is  k,  which  by  superposition  equals 
k'  +  k”.  This  is  the  assumption  upon  which  Eq.  22a  is  based. 

When  Eq.  22a  is  solved  for  y  equals  -i/2,  the  following  expression 
is  obtained  for  the  added  mass  coefficient  of  the  submerged  body, 


k 

s 


d 
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Comparison  of  this  submerged  coefficient  with  that  given  by  infinite  fluid  theory  indicates 
a  divergence  which  increases  with  fineness  ratio  l/d.  In  terms  of  the  present  approxi¬ 
mation  it  will  be  assumed  that  for  y  <>£/ 2  the  infinite  fluid  values  apply  while  for 
y  'K  +  i/2  the  value  of  k  is  zero.  Tne  magnitude  of  the  k  values  given  by  Eq.  22a 
are  adjusted  for  all  values  of  y  by  the  needed  correction.  Thus  the  present  solution, 
as  in  Eq.  22a,  is  utilized  only  in  its  mode  of  variation.  The  variation  in  k  for  motion 
near  the  free  surface  is  finally 

r  -4  i  (-3^2-3(^  +  2(^+|(^,8  +  ...) 


(23) 


The  fully  submerged  value  of  the  added  mass  cocflii  lent  i"*  to  he  determined  from 

potential  fluid  theory  (as  given  in  Fig  Al).  for  t  d  -  4,  k  0  0H2  and  lor  f  tl  - 

k  =  0  029. 

IV  2  Analysis  11 

The  above  analysis  is  subject  to  the  objection  that  the  added* mass  coefficients 
ol  the  elongated  bodies  are  being  determined  as  though  the  bodies  were  spheres  of  the 
same  diameter  as  the  body  cross  set  lion  at  the  free  surface.  An  improved  approximation 
is  to  calculate  the  added  mass  as  that  of  an  ellipsoid  of  revolution  of  diameter  2r  and 
fineness  ratio  (jf/2  -  y).  r  with, of  course,  only  half  of  it  in  the  water.  This  approach 
is  only  valid  for  positive  values  of  y,  as  in  the  previous  case. 

When  the  body  is  less  than  halfway  out  (y  negative!,  the  same  procedure  is 
followed  as  for  the  first  analysis.  Since  it  is  not  possible  to  express  the  relation  for 
the  added  mass  coefficient  in  the  form  of  a  simple  equation  as  with  Eqs.  21a  and  23,  the 
results  are  presented  in  the  Appendix  (Section  VI II  6)  for  two  fineness  ratios.  This 
second  analysis  is  superior  in  requiring  no  correction  to  k  for  agreement  with  the 
fully  submerged  added  mass  coefficient.  The  approx imarion  does  not  seem  completely 
satisfactory,  however,  since  the  added  mass  no  longer  decreases  monotonically  to  zero 
as  the  body  exits.  As  apparent  in  the  tabulation  k,  k^  for  die  finer  body  drops  below 
one  half  before  the  mid  point  of  the  body  is  r "ached,  but  returns  to  this  value  and  tnen 
overshoots  for  awhile.  This  is  the  resulr  of  the  approximation  used,  since  two  segments 
of  the  ellipsoid  placed  end  to  end  do  not  make  another  ellipsoid.  The  single  curve  shown 
in  Fig.  15  was  obtained  by  ignoring  this  hump  and  otherwise  averging  the  results  calcu¬ 
lated  for  the  two  fineness  ratios.  It  is  felt  to  be  an  adequate  approximation. 

An  even  better  approximation  would  be  obtained  by  utilizing  anappt  ox  imate 
formula  for  k,  which  Landweber  and  Wmzer  (2)  found  for  treating  bodies  more  arbi¬ 
trary  than  an  ellipsoid  of  revolution.  This  has  nol  been  tried,  since  the  general  approach 
being  used  is  rather  crude. 

IV.  3.  Analysis  III 

Both  of  the  above  approximations  assume  that  the  variation  in  added  mass 
occurs  only  when  the  body  is  in  contact  with  the  surface.  That  this  may  not  be  so  is 
apparent  from  Eqs.  2  and  2a  and  the  experimental  results  obtained  in  this  study  (1). 

The  added  mass  varies  as  the  submerged  body  approaches  the  free  surface.  At  low 
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Froude  numbers  it  increases,  at  high  Froude  mimi)ers  it  decrease*'  to  as  nun  It  as  is 

nr  SO  percent  of  the  submerged  value 

Since  interest  here  c  enters  in  the  higher  Ftoude  numb*  1  ange  a  ’ational 

improvement  to  the  previous  analyses  appears  possible  by  lssumirg  .*  tie*,  -  ease  in  k 

proportional  to  fy,  f)  *  as  the  body  approaches  the  surface.  Wien  the  hotly  nose  just 

touches  the  surface  it  will  be  assumed  that  k  is  reduced  to  S  S  of  k  .  Actually 

s 

this  reduction  factor  must  depend  upon  »he  Froude  number,  however  the  present 
assumption  appears  to  be  an  improvement  over  that  of  no  variation  in  k  until  the 
surface  is  contacted.  The  resultant  variation  in  k  k^  is  alst  presented  in  Fig.  iS 
and  tabulated  in  the  Appendix. 


4  I . 

V  APPROXIMATE  VKIjOCITY  AMJ  ACCELERATION  EXIT  PROFILES 


O'  i/n-  basis  of  (he  approximate  variation  found  in  (he  previous  chapter  f  > 
the  added  mass  and  relations  for  the  velocity  drag  it  is  possible  to  predict  the  behavior 
of  a  body  during  exit  Of  course,  the  result  will  be  an  oversimplified  picture  ol  whj  t 
happens  since  in  reality  the  exact  mode  of  added-mass-coeff icient  variation  before  and 
during  exit  is  only  approximated.  In  spite  of  this  limitation  it  is  felt  that  these  cal<  ula  - 
lions  are  of  some  value  in  indicating  trends  and  the  relative  importance  of  changes  in 
added  mass,  buoyancy  and  drag  as  the  body  exits.  Furthermore,  they  may  be  of  even 
more  value  when  extended  to  the  angle  exit  case  wh<  re  larger  effects  may  result,  s.  1  n « • 
the  lateral  added-mass  coefficients  involved  are  an  order  of  magnitude  larger 


V .  1 .  Basic  Motion  Equations 

Consider  a  body  of  volume  V,  specific  gravity  S,  and  added  mass 
coefficient  k,  moving  vertically  towards  a  free  surface.  The  forces  acting  on  the 
body  are  buoyance  B,  weight  W  and  resistance  or  drag  D.  Following  Eqs.  11  and 
14,  the  equation  of  motion  for  the  body  is 

W  k 

B  -  W  -  D  -  — (1  +  J)  a  (24) 

R  ^ 

When  the  body  is  at  terminal  velocity  a  =  0  and  B  -  W  -  D  =  In  the  present 

analysis  it  is  assumed  that,  nothing  changes*  until  the  nose  touches  the  free  surface 
After  this  B,  D,  and  k  change.  This  leads  to  a  change  in  velocity  indicated  by 
the  acceleration 


a  _  B  -  W  -  D  =  b  -  S  -  R 
g  WTPHFTS)  ~  $  +  k 


where  b  r  BS/'W  is  a  reduced  buoyancy  force,  varying  from  unity  to  zero  as  the  body 
exits,  and  R  -  DS/W  is  a  reduced  drag  force. 

The  velocity  variation  of  the  body  is  given  as  a  time  integral  of  the  acceler¬ 
ation  rela'ion.  With  the  initial  condition  V  =  at  t  -  0  (well  below  the  surface) 
the  velocity  is 


Except  for  the  Analysis  111  variation  in  k. 
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The  acceleration  term  is  zero  (except  for  k  -  variation  III)  when  the  hotly  is  below  the 
surface  (f  <  •  0)  A  suitable  .space- time  transformation  is  needed  to  effect  this  inte¬ 
gration,  since  b,  k,  and  R  (to  some  extent)  are  known  functions  of  position  of  the  body 
in  the  free  surface  and  not  of  time.  Since  the  change  in  V  is  small  t  is  proportional 
to  y/V ,  whence 

dt  =  l/'V  dy  (27) 

and  the  relation  for  velocity  becomes 

v  -  VT  •  g  /  i  (b  g !■;-■)  dy  (26a) 


The  velocity  V  inside  the  integral  is  either  assumed  constant  or  the  solution  performed 
in  stepwise  fashion 

It  has  been  noted  (cf.  (1).  p.  15)  that  Eq.  23  is  incomplete  when  the  body 
is  near  a  surface  due  to  the  development  of  interference  forces  at  the  surface-  Since 
these  forces  are  usually  not  known,  the  solution  to  the  motion  problem  must  be 
approached  via  energy  considerations.  Lip  to  the  point  where  the  body  breaks  the  free 
surface  it  is  reasonable  to  assume  a  constancy  in  kinetic  energy,  since  the  energy  dis¬ 
sipated  at  the  surface  in  waves  is  probably  small-  Once  the  body  commences  to  actually 
exit  the  change  in  potential  energy  (from  the  motion  of  the  weight)  and  the  work  done  by 
the  buoyancy  and  drag,  which  now  are  variables,  must  be  considered.  With  the  assump¬ 
tion  of  negligible  work  by  the  free  surface  force  (negligible  energy  dissipation  at  free 
surface)  work-energy  considerations  yield  the  following  relation 


(R  -  D  -  W)  Ay  -  A  (kinetic  energy)  =  A 


—  (i  4-  *L) 

|_2g  SW. 


or 

W  k  w 

(B  -  D  -  W)  dy  -  ^(1  +  |)  V  dV  +  2^  dk 

whence 

dV  _  B  -  D  -  W  V  dk 

dy  WV(1  4  k/S)  g  "  2S  (1  +  k/S)  dy 


and  rh..  following  modified  form  of  Eq.  26  results 

~y 


v  -  vT  -  g 


1  ,b  -  S  -  R,  .  1 

v  (TTF)  (|y  '  1 


V 


(s+k) 


dk 


(28) 


It  appears  that  Eq  26a  is  but  a  first  approximation,  depending  upon  the  magnitude  of  the 
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.second  term  in  Eq  28  li  spec  non  suggests  that  since  ns  relative  significance  v.uies  as 

2 

kV\  this  term  may  tlius  he  (|uite  important  for  rathei  blunt  bodies  jf  appreciable  Fioude 
numbers. 


V  2 .  Relation  for  Reduced  Buoyancy  Force 

Before  calculating  the  acceleration  and  velocity  variations  of  the  rn'ing 
bodies  it  is  desirable  to  organize  information  on  buoyancy  and  drag  variations  in  suitable 
fashion.  The  relations,  b  =  b  (y/jf)  and  R  ~  R(y/f)  are  needed.  The  third  variable 
in  Eq.  26  or  28  is  k  “  k(y  ()  and  this  is  available  from  the  preceedmg  section. 

The  relation  between  b  =  BS/W  and  the  position  of  the  body  in  the  inter¬ 
face  is  simply  one  of  geometry.  When  the  body  is  partly  exited,  its  buoyancy  is  s.mply 
that  of  the  residual  submerged  volume  ¥’  times  the  liquid  specific  weight.  Since 
W  =  Sy  V,  where  y  is  the  specific  weight  of  water  and  V-  the  rc  fa*  volume  of  the 


body 


.  RS  VyS  ¥' 
b  ‘  W  '  5TV  '  Y 


(29) 


From  the  geometry  of  the  el lip.se> id,  with  the  aid  of  calculus,  the  volume  ratio  is  given 
by  the  relation 


b 


¥_ 

¥ 


(29a) 


As  In  previous  sections  y  is  the  distance  of  the  center  of  the  ellipsoid  from  the  free 

1  V  A 

surface.  Equation  2Ha  app  :es  only  in  the  range  ~  J  <  J  <  +  ~2  w^f'>n,  e  ^  vanes 

from  unity  to  zero  as  shown  in  Fig.  16  (cf.  page  39). 


V.3.  Reduced  Resistance  Force 

The  drag  depends  upon  the  wetted  surface  of  the  body  and  its  velocity 
through  the  water.  In  contrast  to  the  calculations  of  Chapter  II,  the  wetted  or  surface 
area  is  now  a  variable  Moreover,  it  has  been  noted  in  the  experiments  that  when  the 
rising  body  is  partly  out  of  the  wafer  an  appreciable  amount  of  water  clings  to  its  sur¬ 
face.  Considerations  of  this  later  effect  will  be  deferred  to  the  next  chapter  and  for  the 
present  It  is  assumed  that  the  body  above  the  interface  is  dry.  Since  air  is  much  less 
dense  than  water  its  drag  will  be  neglected.  It  is  also  assumed  that  the  velocity  varia¬ 
tion  of  the  body  during  exit  is  relatively  small,  so  that  the  drag  relation  R  -  DS/W 
=  D/yV  is  a  unique  function  of  y/f. 
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Determination  of  the  reduced  resista'i  <•  ri  Lit  ion  .m  o1'. «  evaluation  of 
the  manner  in  which  the  skin  friction  and  pressure  portion***  >,i  the  d'ag  vary  along  the 
length  of  the  body  IJp  to  the  ;>oint  of  boundary- layer  separation  the  pi sse*e  distribu¬ 
tion  along  the  body  must  be  i  lose  to  that  obtainable  fron  ideal  -  fluid  theory  Pie  jtoitit 
of  separation  and  sk  in-friction  distribution  are  determined  by  the  na'\re  ,,f  -|)L.  }>oundary- 
layer  development  along  llie  body  For  the  laminar  boundary  layer  an  u.  the  separation 
point,  thts  has  been  determined  in  general  form  by  PretscL  (17;  based  <  r  the  ideal- 
fluid  velocity  distribution  along  the  surface.  The  calculation  c!  the  vacation  along  die 
surface  of  the  drag  contribution  of  die  rwo  factors  is  thus  relatively  straight  forward, 
up  to  the  separation  point. 

Figure  17  depicts  the  pressure  distribution  along  the  body  of  fineness 
ratio  4,  in  terms  of  the  pressure  coefficient 


where 
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p  -  pressure  along  die  surface 

p^  =  refercnce  pressure  measured  upstream 

IJ  -  velocity  along  body  surface  (just  oe'side  th<  boundary- layer) 

(=V)  -  reference  velocity  (relative  velooi'y  between  the  body  and 
undisturbed  fluid) 
p  -  fluid  density 


The  ideal-fluid  pressure  distribution  is  symmetrical  fore  end  aft.  Jn  the  real-fluid 
case,  with  separation,  the  pressures  near  the  tail  are  somewhat  less  as  will  be 
discussed  later-  The  variation  in  skin  friction,  as  taken  from  Prefsch  s  work.,  is  indi¬ 
cated  by  the  parameter 

t~-c./k~  (.11) 


where 


*  Calculation  of  the  pressure  drag  variation  along  the  ellipsoid,  in  similar  fashion 
to  the  present,  was  given  by  A.  F  Zahm  in  192b  (NACA  Report  253)  both  from  the 
ideal-fluid  pressure  distribution  and  the  pressures  measured  by  Jones  (14). 
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C 

r 


skin  ft  iction  coefficient 


o 

r  =  local  wall  skin  friction 
w 

=  Uq  t  v  -  Reynolds  number  of  body  in  fluid 

This  parameter  T  is  independent  of  Reynolds  number,  except  in  regard  to  occurrences 
after  separation  For  these  calculations  z  is  taken  as  the  distance  along  the  major 
axis  from  the  body  nose  and  r  as  the  lateral  coordinate  of  the  surface  !n  terms  of 
ds,  the  element  of  length  along  the  surface  and  the  tangent  angle,  a  =  sin  dz/ds 
=  cos  1  dr/ds,  the  surface  forces  at  any  position  along  the  major  axis  are 


dD 

s 

dD 

P 

where  the  subscri 
fions  of  the  drag. 


D  = 


=  2  tr  r  r  cos  a  ds  =  2  tt  t  rdz 

w  w 

=  2  tr  r  p  sin  a  ds  =  2  tt  p  r  dr 

pts  s  and  p  indicate  the  skin  friction  and  normal  pressure  por- 
Integration  of  these  over  the  body  will  yield  the  body  drag 


Reducing  this  relation  to  dimensionless  form,  through  division  by  the  dynamic  pressure, 
yields  the  relation 


'D 


2  D 

P  U  2  A, 
o  J- 


(32) 


r  ,  ,r, 


cPid<? 


where  is  the  normal  area  of  the  body  (projected  on  plane  perpendicular  to  Uj) 

and  f  is  the  length  diameter  ratio  of  the  body  (4  in  the  case  calculated) 

For  the  water-exit  calculations  it  is  necessary  to  evaluate  the  integrals  in 
Ec].  32  at  various  positions  along  the  body  from  0  to  z/f  -  1  Thus  for  z/£  =  0  50 
this  will  indicate  the  drag  contribution  of  the  front  half  of  the  body.  Dividing  the  results 
of  such  partial  solutions  of  Eq.  32  by  the  result,  yields  die  relative  contribution 

to  Uie  drag  of  various  portions  of  the  body,  forward  of  any  z/jf  station  The  result, 
calculated  with  the  aid  of  the  c:  and  T  values  of  Fig.  17,  is  shown  in  Fig  IK. 


Cumulative  Drag  Contribution  A,  D 
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From  (Jus  figure  it  appears  that  the  rotal  diag  contributions  of  the  [xirtion  of  the.-  Ixnly 
forward  of  a  given  station  may  be  as  low  as  45  percent  ot  as  high  at  22H  peuent  of  the 
resultant  body  drag  IXj.  This  variation  is  primarily  due  to  the  pressure  For  use  in 
the  water-exit  calculation,  it  will  be  assumed  that,  as  the  body  exits  the  drag  lontn- 
butions  of  the  still-wetted  portions  remain  the  same  .is  they  were  when  the  Ixuly  was 
completely  submerged.  The  residual  drag  contributions  of  the  wetted  portions  are 
thus  found  by  subtracting  the  contribution  indicated  by  Fig.  »8.  at  the  appropriate 
z,  l  station,  from  unity.  The  relation  for  R  will  be  presented  in  )i(\  33. 

Before  proceeding  with  the  water-exit  calculations  a  few  remarks  about 
the  post-separation  behavior  of  the  flow  along  the  body  surface  are  in  order.  The 
boundary- layer  analysis  gives  no  clue  as  to  occi  rences  downstream  from  the  station 
z/ ft  -  0.838  where  separation  is  indicated.  The  flow  reversal  after  separation  leads 
to  negative  values  of  T,  but  these  have  little  influence  on  die  total  drag  --  most  of 
the  skin-friction  drag  is  produced  on  the  front  of  the  body.  The  assumed  T  variation 
near  the  tail  of  the  body,  as  depicted  in  Fig.  11,  is  felt  to  be  reasonable.  The  flow 
behavior  after  separation  has  a  much  more  significant  effect  on  the  pressure  drag. 

The  values  for  Cp  used  in  the  calculations  near  the  tail  and  shown  in  Fig.  17 
represent  a  considerable  regain  in  pressure  following  separation.  Since  little  informa¬ 
tion  on  the  flow  behavior  and  resultant  pressure  regain  after  laminar  boundary- layer 
separation  was  found  in  the  literature,  the  regain  was  adjusted*  to  yield  a  reasonable 
drag  coefficient.  From  Fig.  4  for  an  ft/ d  =  4  body  without  fins  and  no  transition  to 
turbulence  effects,  the  drag  coefficient  af  IR^  =  10  is  -  0.0032  In  terms  of  the 
fineness  ratio  of  the  body  Cp  is  0.  0412  (i.  e  12.87  ) 

The  drag  contributions  of  the  fins  used  on  the  test  bodies  must  still  be 
included.  Since  these  are  simple  flat  plates  the  calculation  is  straight  forward  In 
terms  of  the  ftn  proportions  for  the  2- in  £/d  -  4  body  the  result  is  depicted  in  Fig 
18  as  the  dashed  line.  The  variation  in  A"D^  >s  thus  from  zero  to  1. 43  to  0  33  at 
about  the  midpoint  on  the  body,  to  1  92  and  back  to  1.00  at  the  tail. 


*  A  regain  in  pressure  to  a  C  of  aboui  0.45  at  the  rear  stagnation  point  was  required. 
The  effect  of  the  pressure  ^  regain  on  drag  is  large.  Thus  the  first  estimate  of  a 
regain  to  0.  10  yielded  a  Cp  value  of  0  0816,  nearly  twice  the  value  deemed  appro¬ 
priate.  The  observations  J  of  [ones  (14)  on  an  f  -  4  ellipsoid  were  particularly 
suggestive  in  this  regard  He  found  Dp  to  be  40  percent  of  1)  versus  the  present 
result  of  55  perc  ent. 
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where  IXj.  is  evaluated  from  the  terminal  veloe 1 1 \  (below  surface  influence  region) 
relation  of  Lq  8.  It  should  be  noted  again  that  changes  in  the  velocity  of  the  body  arc- 
assumed  to  be  small  enough  to  have  negligible  effect  on  R.  With  13.^.  varying  with 
velocity  to  some  exponent  usually  less  than  2  and  in  terms  of  the  ex  pi"  imental  velocity 
changes  of  the  order  of  5  percent  or  less,  the  error  introduced  by  tins  approximation 
is  not  large. 


V.  4.  Velocity  and  Acceleration  Profile  Calculations 

Solution  of  Kq  28  via  stepwise  integration  with  the  aid  of  the  relations 
developed  for  b,  R  and  k  has  been  carried  out  for  a  number  of  situations  Calcula¬ 
tions  have  been  limited  to  the  f  -  l,  d  -  4  ellipsoids  for  which  fc  f  Fig.  At) 

k  -  0.  08  1  f>. 
s 

Figure  19  presents  some  of  the  results  obtained  for  a  body  of  specific 
gravity  S  =  0.  25  at  a  Froude  number  F  -  8.  These  are  presented  in  the  same 
fashion  (cf  Fig.  3)  as  used  for  (lie  experimental  results  (1),  i  e.  as  the  velocity  ratio 
V/VT  and  acceleration  a  versus  s/9 }  the  relative  position  of  the  center  of  the 
body.  Portions  of  the  acceleration  profiles  are  shown  dashed  in  regions  where  the 
stepwise  integrations  were  not  detailed  enough  (the  increment  of  calculation  ^y  H 
was  too  large)  to  clearly  delineate  the  mode  of  variation  in  a.  Those  curves  repre¬ 
sent  smoothed  lines  through  the  calculated  values.  For  the  body  center  below  the  sur¬ 
face  the  analytical  results  bear  considerable  resemblance  to  the  experimental  Above 
the  surface  the  second  hump  in  the  velocity  and  acceleration  profiles  does  not  appear 
experimentally .  One  possible  factor  (entra  ined  water)}whieh  may  adjust  for  this 
effect,  is  discussed  in  the  following  chapter  The  present  discussion  will  therefore  be 
limited  to  fhe  below  surface  results.  Before  leaving  the  above  surface  predictions;  it 
is  interesting  to  note  that  a  peak  velocity  ratio  of  about  i.  19  is  predicted,  appreciably 
less  chan  the  value  1 .  .33  predicted  by  the  elementary  momentum  analysis  (velocity  jump) 
as  indicated  by  Fq.  1.  The  below  surface  peak  indicated  in  Fig.  19  is  1.11  to  1.  14, 
less  than  the  above  surface  peak. 


Fig-  19.  Calculated  Velocity  and  Acceleration  Profiles  for  Different  Variations  in 
Added  Mass. 
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A  major  objective  of  the  analysis  summarized  here  was  to  evaluate  the 
relative  significance  of  the  various  factors  entering  into  the  velocity  and  acceleration 
profiles.  One  question,  which  arises  when  more  exact  analysis  of  the  water-exit  pro¬ 
blem  is  contemplated,  is  whether  a  refined  analysis  of  the  variation  in  added-mass  co¬ 
efficient  is  worthwhile.  It  is  possible  that  the  effects  of  buoyancy  and  resistance  change  s 
on  the  behavior  of  the  exiting  body  may  far  outweigh  the  effects  of  added-mass  changes. 

A  partial  answer  to  this  question  is  evident  from  the  results  of  the  calculations  pre¬ 
sented  in  Fig.  19.  Here  four  different  variations  in  added-mass  coefficient  have  been 
employed  in  the  calculations,  namely,  that  marked  0,  with  k  constant  (k  =  k ^ ) , 
and  the  I,  II,  and  III  variations  of  Chapter  IV.  The  results  based  on  Analysis  III  for 
k  approach  closest  to  the  experimental  results  where  velocity  and  accelerations  in¬ 
creases  are  observed  to  commence  before  the  body  nose  reaches  the  surface  (y,  ( 

=0.5).  The  profiles  in  this  below-surface-contact  region  result  entirely  from  the 
second  term  in  Eq.  28.  The  holding  of  k  invarient  throughout  the  water-exit  period 
(curve  0)  considerably  reduces  the  velocity  changes  and  accelerations.  Thus  the  below- 
surface  peak  is  only  1.042  in  this  case,  rather  close  to  the  experimentally  observed 
value  Comparison  of  the  several  k  relations  used  for  the  calculations  presented  in 
Fig.  19,  thus  indicates  that  the  variation  in  k  has  a  significant  effect  but  that  other 
factors  are  far  from  insignificant.  Unless  much  more  sudden  changes  in  k  appear 
probable,  the  exact  details  of  the  variation  are  not  important.  Thus  curves  I  and  H 
differ  little. 

The  double-humped  feature  of  the  calculated  profiles  results  from  the 
variation  in  the  resistance  factor  R  (cf.  Eq.  33)  and  A/D^  (cf.  Fig.  18)  with 
y/f.  Calculation  of  the  profiles  which  would  occur  with  R  invarient  are  presented 
in  Fig.  20.  The  velocity  and  acceleration  are  seen  to  peak  quickly  and  then  to  drop 
rapidly.  The  resemblance  to  the  experimental  profiles  (no  second  hump,  etc.)  is 
closer  except  that  the  calculated  decelerations  become  much  larger  than  the  observed 
ones. 

The  influence  of  body  specific  gravity  S  and  Froude  number  IF  on 
the  water  exit  phenomena  have  been  indicated  in  Chapter  1.  Thus  Eq.  1  suggests  that 
the  velocity  increase  should  vary  inversely  with  S.  The  experiments  tend  to  verify 
this  as  well  as  to  indicate  a  Froude  number  effect  (cf.  Fig.  4).  Reference  (1)  shows  a 
rapid  increase  in  Froude  number.  Calculations  of  the  profiles  have  also  been  carried 
out  for  two  Froude  numbers  and  two  specific  gravities.  The  results  are  presented  in 
Fig.  21  which  appears  to  verify  the  expected  effect  of  S  on  the  velocity  increase  as 
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ci  Acceleration  Profiles. 


well  us  to  indicate  un  effect  of  S  on  the  accelera t ions .  The  effect  i.s  not  well 

correlated,  however.  An  increase  in  acceleration  with  IF  is  found  sum  what  as  m 
the  experiments  hut  (lie  variation  in  relative  velocity  increase  with  IF  shown  in  Fig. 
21  is  not  found  experimentally. 

'Hie  assumption  of  a  lack  of  Froude  number  effect  on  k  (cf  Section 
IV  3)  presumably  enters  here.  These  calculations  have  been  limited  to  the  f  -  f/d 
-  4  body  shape.  'I3ie  major  effect  of  fineness  ratio  f  is  a  change  in  the  submerged 
added  mass  coefficient  k 
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VI.  CONSIDERATIONS  OF  ENTRAINED  WATER 

From  ihe  first  tests  (in  November  1956),  the  experimental  phases  of  this 
project  have  ind  cated  that  significant  amounts  of  water  were  entrained  by  the  exiting 
bodies.  Figure  2  indicates  this  clearly  Crude  measurements  of  this  quantity  have 
been  reported  (1)  and  the  volume  was  shown  to  be  an  appreciable  fraction  of  the  volume 
of  the  body  --of  the  same  order  of  magnitude  as  the  added-mass  coefficient.  It  appears 
probable  that  this  entrained  water  is  associated  with  the  boundary  layer  on  the  body. 
However,  die  displacement  volume  of  this  layer  was  calculated  and  found  to  represent 
only  about  55  percent  of  the  measured  entrained  water  (1).  These  calculations  are  pre 
sented  m  more  detail  in  the  following  section  after  which  they  are  extended  to  other 
scales  of  motion  in  an  attempt  to  assess  the  possibilities  of  Reynolds-number  scaling 
on  entrained  water.  In  the  previous  chapter  mention  was  made  of  the  possible  effect 
of  this  entrained  water  on  the  water-exit  calculations.  This  factor  is  introduced  into 
the  motion  analysis  in  the  final  section  of  this  chapter. 

VI.  1  Calculation  of  Displacement  Volume 

The  growth  in  displacement  thickness  6*  along  the  body  may  be  c.alcu- 
lated  from  the  work  of  Pretsch  (17).  He  presents  curves  for  various  ellipsoids  depicting 
the  variation  in  VlR ^  6*/|  along  the  body  surface.  One  is  shown  in  Fig.  17,  for  the 
f  =  Sj d  =  4  body  shape,  versus  the  coordinate  z  along  the  body.  The  displacement 
volume  •¥■*  of  the  boundary  layer  is  found  from  6*  by  the  following  integral 


In  terms  of  a  dimensionless!  reduced  form,  appropriate  to  laminar-boundary- layer 
analyses,  this  may  be  written 

A 

v'sr,  *  2  y  Jr]  \  d  \  (.34) 

where,  as  before,  i  is  th«.  body  length,  r  the  radius  of  its  cross  section  at  any 
axial  distance  z  from  the  nose.  If  the  boundary  layer  dtd  not  separate  (cf.  Section 
V.  3),  Eq.  34  could  be  integrated  to  the  tail  of  the  body  (z/f  -  1)  to  yield  the  displace¬ 
ment-volume  relation.  However,  separation  does  occur,  at  a  z/f  value  of  about  0  H  ^ 


(cf.  Fig  17)  and  in  (he  wake  beyond  this  point  Jli  ^  f>*  4  values  are  n<  i  available 

from  boundary- layer  analysis. 

The  boundary- layer  displacement  volume  up  to  the  serration  point  is 
easily  calculated  from  Pretsrh's  result  introduced  into  Eq.  34.  The  variation  in  the 
integrand  with  z /  f  is  nearly  linear,  so  the  calculation  is  simple.  The  boundary- 
layer  displacement  volume  for  a  laminar  boundary  layer  on  an  ellipsoidal  body  of 
f  =  4,  to  a  point  (z/j ?  =  0.  82)  just  before  separation,  is  found  to  be 


=  0.510//!^ 


or  in  terms  of  the  body  volume  V--  0.0327  4  (cf.  Fig.  Al), 


V-*\ 

£  =  15-6  (1R  ) 


1/2 


(3S) 


where  the  subscript  4  is  used  to  indicate  the  laminar  layer  and  b  the  body.  Thus, 

at  a  Reynolds  numoer  of  10^  the  Laminar-boundary- Layer  displacement  volume  is  1.  56 

4 

percent  of  the  body  volume  and  at  10  it  is  15.  6  percent. 

The  wake  contribution  V*  to  the  entrained  displacement,  volume  can  only 

^  -f 

be  estimated.  The  behavior  of  the  boundary  layer  after  separation  is  not  well  known, 
as  has  been  noted  previously  (cf.  Section  V.  3).  This  flow  is  often  termed  the  laminar 
burble,  but  it  rapidly  becomes  the  wake.  At  Reynolds  numbers  of  interest  here  (of  the 
order  of  10^)  this  separated  flow  will  not  remain  laminar  for  long.  After  transition  to 
turbulent  flow  such  separated  flows  have  a  tendency  to  reattach  (18).  In  the  present 
case  the  separation  point  is  near  the  tail  where  the  body  is  so  blunt  that  it.  has  been 
assumed  there  is  no  reattachment.  Further  consideration  of  the  meager  information 
(18,  19)  on  the  laminar  burble  suggeststhat  the  flow  leaves  the  surface  along  a  tangent. 
The  wake  displacement-volume  increment  was  therefore  calculated  as  a  conical  frustrum 
tangent  to  the  displacement- thickness  contour  just  before  separation  Th>s  volume  was, 
of  course,  reduced  by  the  volume  of  that  portion  of  the  body  contained  within.  The  cal¬ 
culation  was  terminated  at  the  tail.  This  calculation  was  made  separately  for  each 
Reynolds  number  of  interest,  since  the  relative  boundary-layer  thickness  varies  with 

Reynolds  number.  However,  the  wake  relative- volume  contribution  -V*  /V-  was 

6  w 
essentially  a  constant.  At  10  the  wake  contribution  was  0.  035  or  some  70  percent  of 

the  total. 


+  Technically  after  separation  the  boundary  layer  indy  r.:: 


|nnije'  be  so  termed. 


The  (mill  Hem  to  he  included  is  die  viscous  displacement  volume  genet  a  ted 

by  die  (our  fins.  Hie  displacement  thickness  growdi  on  these  is  easily  computed  from 

flat-plau relations  At  a  Reynolds  number  of  10^  the  (in  contribution  V*/V  is  about 

0.002  or  4  percent  In  general  the  fin  contt  ibution  is  2  /* F  .  •  At  the  same  Reynolds 

5  * 

number  (actually  9.  7  x  10  )  the  total  boundary- layer  displacement  volume  was  foend  to  be 


The  measured  entramed-water  volume  for  the  2-in.  body  was  1.7  times  this.  In  the 
case  of  the  1-in  body  (at  2.  9  x  10^),  the  measured  water  was  1.85  times  the  calculated 
displacement  volume  (1).  Thus,  the  correlation  between  bounda ry- layer  displacement 
volume  and  measured  entrained  watei  is  only  fair. 

'rile  divergence  noted  between  measured  entrained  water  and  calculated 
boundary- layer  displacement  volume  may  be  due  to  any  one  of  several  factors.  Thus 
the  experimental  n  easurements  were  crude.  They  may  have  included  a  considerable 
amount  of  water  in  ihe  wake  behind  the  body  and  this  latter  was  not  calculated.  Or. 
surface  tension  or  other  fluid  properties  and  phenomena  may  affect  the  entrained  water. 

VI.  2.  Displacement  Volume  for  Other  Scales 

In  spite  of  the  uncertainties  and  apparent  inaccuracies  in  the  correlation 
between  entrained  water  and  displacement  volume  it  seems  useful  to  estimate  the 
variation  of  this  volume  with  the  scale  of  test,  on  the  assumption  that  this  indicates  the 
manner  with  which  the  entrained  water  will  change  with  scale.  The  experimental  pro¬ 
gram  of  water  exit  studies  (1)  has  involved  bodies  of  diameters  from  1  to  4  inches. 

This  is  a  considerable  range  of  sizes  and  of  course  there  is  always  the  question  of  ex¬ 
trapolating  experimental  results  to  possibly- larger  prototype  sizes. 

The  calculations  summarized  in  the  preceding  section  were  performed  at 
Reynolds  numbers  corresponding  to  test  results  with  the  1-in.  and  2-in.  bodies.  Thus 
at  2.  9  x  10S  die  calculated  ¥-*/¥  value  was  0.  068  and  at  9.  7  x  1 it  was  0.  053.  The 
measured  entrained  water  similarly  increased  with  decrease  in  scale.  The  question  is, 
what  are  the  values  for  the  4- in.  body  at  3.  5  x  10^  or  some  other  even  larger  body? 

If  boundary- layer  separation  and  transition  at  the  higher  Reynolds  numbers  do  not  enter 
to  complicate  the  flow  pattern,  a  relation  similar  to  Eq.  34  would  provide  a  simple 
direct  answer.  The  relative  displacement  volume  would  vary  inversely  as  the  square 
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Fig.  22.  Variation  in  Relative  Displacement  Volume. 


root  of  the  Reynolds  number*.  The  wake  contribution  is  not  simply  calculated  once 

and  for  all,  even  without  transition.  According  to  Schlichting,  although  the  point  of 

laminar-boundary-layer  separation  on  a  given  body  is  Invariant  with  Rt7nolds  number, 

the  angie  with  which  the  flow  leaves  the  surface  is  IR  dependent.  There  appears  to 

be  little  information  on  just  how  this  angle  varies,  however.  The  wake  contribution 

has  thus  been  evaluated  as  in  the  previously  mentioned  calculations.  The  relative 

displacement  volume  estimate,  for  the  laminar  boundary- layer  is  shown  in  Fig-  22  as 

4 

a  function  of  Reynolds  number.  It  is  seen  to  decrease  from  22  percent  at  10  to  about 


For  Reynolds  numbers  above  4  x  10  since  below  this  value  on  a  flat  plate  the 
Blasius  solution  no  longer  is  correct. 


1 


* 


ve 

•1  percent  al  1U 

At  the  higher  Reynolds  ininihfi!)  (alx>ut  1  a mJ  aixjve)  i(  becomes  increas¬ 
ingly  piohablc  l ha t  the  Ixiundary  layer  on  the  !>jdy  will  become  turbulent  before  the 
laminar  separation  jxnnt  is  reached.  This  expectation  is  indicated  by  die  expected  drag 
coefficient  curves  of  l;ig.  3.  with  some  experimental  verification  shown  in  Fig.  9.  The 
turbulent  boundary  layer  being  more  resistant  to  separation,  transition  will  thus  lead 
to  a  lower  wake  contribution  to  the  bounda ry- layer  displacement  volume.  5ince  this 
term  represents  70  or  more  percent  of  die  total  at  die  higher  Reynol  'i  numbers,  the 
effect  of  transition  in  reducing  the  displacement  volume  is  considerable. 

On  die  basis  of  a  transition  (to  turbulent  boundary  layer)  Reynolds  number 
of  10^  the  displacement  volume  has  been  estimated  at  body  Reynolds  numbers  of  2  x  10^ 
and  4  x  10^.  The  calculation  is  b.e.cJ  on  use  of  the  Karman  momentum- integral  equation 
in  step-wise  calculation  from  die  transition  point  The  method  used  is  essentially  that 
given  by  Granville  (12)  except  that  corrections  were  introduced  into  the  boundary- Layer 
shape  parameter  and  skin  friction  following  the  work  of  Ross  (20).  Separation  was  in¬ 
dicated  at  a  z  !  value  of  about  0.  95,  in  contrast  to  the  laminar  separation  point  of 
0.  838.  This  led  to  a  considerable  reduction  in  die  wake  contribution,  while  the 

boundary-  layer  portion  is  higher,  as  may  be  seen  in  Fig-  22.  As  suggested  in  the 

t) 

figure,  at  Reynolds  number  above  4  x  10  ,  the  displacement  volume  must  reach  a  peak 
value  and  then  commence  to  decrease,  after  the  boundary  layer  on  the  wnole  length  of 
the  body  becomes  turbulent.  The  sudden  drop  in  -V-*/V  at  about  1.  2  x  10^  occurs 
when  the  transition  point  changes  from  being  in  the  separated  flow  to  being  on  the  body 
and  the  separation  is  then  delayed. 

Referring  these  predictions  of  displacement  volume  to  the  actual  water- 
exit  test  conditions,  it  is  notable  that  for  the  4-in.  body  tested  (at  IR  ^  =  3.4b  x  10^) 

¥■•/¥■  has  a  calculated  value  of  2.  7  percent  versus  5.  2  percent  for  the  2- in.  body  (at 
10^)  and  6.7  percent  for  the  1-in.  body  (at  3  x  10^).  Whether  this  properly  represents 
or  indicates  the  entrained  water  or  not,  it  is  obvious  that  a  considerable  scale  effect 
does  occur.  This  is  part icula rly  so  at  sizes  above  the  2-in.  body  in  association  with 
boundary- layer  transition  effects. 

VI.  3.  Hntrained  Water  Effect  on  Exit  Profiles 

With  an  amount  of  water  exiting  with  the  body  of  about  the  same  value  as 
the  added-mass  coefficient  a  cons idedable  modification  of  the  exit  behavior  is  expected. 
The  question  of  whether  this  viscous  displacement  volume  simply  increases  the  added 


60 


mass  coefficient  k  ,  involves  consideration  of  unsteady  bounda ry  -  layer  behavior  A( 
exit  however,  the  velocity  changes  are  small  and  a  considerable  amount  of  water  is 
observed  to  exit  with  the  body.  It  will  here  be  assumed  that  this  displacement  volume 
actually  clings  to  the  body  during  the  entire  exit  phase.  It  has  been  observed  in  the 
motion  pictures  that,  when  the  body  is  most  of  the  way  out  of  the  water  some  water 
commences  to  roll  down  and  off  the  body.  This  suggests  that  the  boundary- layer  dis¬ 
placement  water  does  cling  to  the  body  during  most  of  the  exit  phase.  An  attempt  to 
evaluate  the  effect  of  this  water  on  the  calculated  velocity  and  acceleration  profiles 
(cf.  Section  V.  4)  has  been  made  and  will  be  presented  in  the  following  paragraphs. 

It  is  assumed  that  the  boundary- layer  displacement  volume  clings  to  the 
body  and  exits  with  it.  This  Implies  that,  although  the  buoyancy  change  is  the  same 
during  exit,  the  mass  of  body  plus  water  exiting  is  larger;  also  while  the  boundary  Layer 
is  thin  enough  for  the  variations  in  pressure-drag  contributions  to  be  the  same,  the 
frictional-drag  contribution  will  not  vary  as  the  body  exits.  Thus  considerable  change 
is  introduced  into  the  several  terms  in  the  equations  leading  up  to  Eq  28,  particularly 
R  and  W-  Splitting  R  into  pressure  and  skin  friction  parts  R^  and  Rg  with 
the  latter  invarient  with  y / i  and  introducing  the  relative  displacement  volume 
v*  =  ¥*/V  into  appropriate  places  in  the  several  terms,  Eq  28  now  appears  as 


b  -  R  -  R  -  S  -  v4 
/I  +  vVSKS'Tl)  ' 


dk 


(28a) 


The  numerator  of  the  bracketed  term  in  the  first  integrand  contains  the  extra  term 
v*  and  only  b  -  Rp  is  a  function  of  body  position  at  the  free  surface.  The  second 
integral  is  not  changed  at  all.  The  b  -  Rp  term  is  available  from  Fig.  16  for  b 
and  Fig.  18  for  Rp  This  latter  term  is  42  percent  of  the  total  drag  or 


R  =  0.42  (1  -  A  /D  )  (1  -  S) 
P  P  T 


(33a) 


where  Ap/D^  is  found  from  the  pressure  curve  of  Fig.  18  multiplied  by  the  factor 
0.  763  =  l/l.  31,  a  correction  for  the  fact  that  this  pressure  curve  is  presented  as  a 
fraction  of  the  bare  body  drag  and  the  fins  added  31  percent  to  the  bare  body  drag.  The 
R  term,  representing  58  percent  of  D_  is  simply  0.581  (I  -  S),  a  constant. 

S  1 

Velocity  and  acceleration  profiles  have  been  calculated  on  this  new  basis  in 
the  same  stepwise  fashion  as  used  previously  (cf.  Section  V.  4).  Results  are  presented 
in  Fig.  23,  in  comparison  with  the  previous  results  for  one  Froude  number  and  body 
density.  The  below-surface-contact  behavior  is  naturally  the  same  as  found  previously 


.  23.  Calculated  Exit  Profiles  With  and  Without  Estimated  Wated  Entrainment. 


since  this  results  only  from  the  second  integral  of  Hq  28  which  is  not  (.Linked  An 
appreciable  reduction  in  the  velocity  jump  and  peak  acceleration  is  seen  to  result 
from  the  water  entrainment  effect.  Furthermore  the  second  peak  in  veil  uty  and 
acceleration  has  also  disappeared.  The  calculated  profiles  thus  more  closely  resemble 
those  observed  (cf.  Fig.  3).  The  discrepancy  is  now  mere  one  of  degree  A  greater 
entrained  water  effect  could  lower  the  calculated  velocity  profiles  m  the  regions 
-0.2  <  y,  (  until  rather  good  agreement  with  experiment  is  found.  The  question  of 
the  absolute  accuracy  of  the  entrained  water  calculations  as  the  boundary  layer  dis¬ 
placement  volume,  in  light  of  the  experimental  measurements,  has  already  been 
discussed  in  Section  VI.  1. 

Further  refinement  in  the  calculation  of  the  profiles  seems  inappropriate 
at  this  time.  These  last  calculations  have  definitely  shown  a  significant  and  appreci¬ 
able  effect  of  the  water  entrainment  on  the  profiles.  If  this  entrained  volume  varies 
with  the  scale  of  test  (as  indicated  by  Reynolds  number)  in  the  fashion  suggested  by 
the  volume  calculations  summarized  in  Fig.  22,  then  Reynolds  number  scaling  also 
influences  water  exit. 
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as  a  suhsc  -  ipi  indiLd r«  s  ..  i r 

r  A  t  d  -  body  su^jh.'  ,i  i  cu  •- elated  in  normal  dimensions  of  body 

we'fed  (.  -  surfai  e  a  ca  c  1  body,  ft" 

> 

fin  area  ft" 

nor  mal  or  projec  led  a : ea  ft " 

=  BS,  W  -  reduced  buoyancy  force 
buoyancy  force  of  submerged  body,  lbs 
drag  coefficient  based  or  A  .  cf.  Eq.  1 
flat-plate  friction  coefficiert 
pressure  coefficient,  cf.  Eq.  30 
local  skin-friction  coefficient 
diameter  of  body,  inches  or  feet 
body  drag,  lbs. 

body  drag  in  infinite  fluid  at  terminal  velocity 
body  drag  due  'o  skin  friction 
body  drag  due  to  pressure 
=  f/d  -  fineness  ratio  of  body 

-  V/  /gd  -  Froudc*  number 
gravita'ic  nal  acceleration,  32.  2  b  s^ 
added-mass  coefficient 

added-mass  coefficient  for  body  submerged  in  infinite  fluid,  Chapter 
IV,  V 

coefficient  in  power-law  relation  for  drag,  cf.  Eq.  12 
length  of  body  ft. 
mass  cf  body  slugs 
added  mass  of  body,  slugs 

exponent  in  power- law  relations  for  drag,  cf.  Eq.  12 

pressure  along  body  surface  psf 

refe’om  c  or  ambient  pressure  ahead  of  body,  psf 

tad  us  of  body  fperpendicula 1  to  major  axis)  at  any  point,  ft 

-  1  )S  /  W  -  ’educed  drag  factor.  c.t.  Chapter  V 


-  V  tj  v  -  Reynolds  number  based  on  body  length 

specific  gravity  of  body  =  ra»i<  of  body  weight  to  that  of  equal 
volume  of  water 

time,  seconds 

-  C  </  IR  .  -  reduced  lamina'  friction  coefficient 

T  ft 

velocity  along  body  surface  just  outside  boundary  layer  (potential 
velocity),  fps 

-  V  -  relative  velocity  of  body  through  fluid,  fps 

3 

-  V/  £  -  body  volume  related  to  length  cubed 
=  V*/V  -  relative  displacement  volume 
velocity  of  body,  fps 

terminal  velocity  of  body,  fps 

3 

volume,  usually  of  or  pertaining  to  body,  ft 

3 

boundary- layer  (etc. )  displacement  volume  on  body  ft 
as  a  subscript,  indicates  water  except  in  Section  v  3 

-  ySV  -  weight  of  body ,  lbs 

V/V^  -  relative  velocity  in  Chapter  IM 

position  of  center  of  body  relative  to  free  surface,  ft 
(negative  when  below  surface,  positive  above) 

/  2 

yg/a/JV^  -  rise  distance  parameter  in  Chapter  III 
position  along  major  axis  of  body,  measured  from  nose,  ft 
=  1  +  k/S 

surface  tangent  angle 
=  S/(1-S) 

fluid  specific  weight.  lbs/fr'^ 

boundary- layer  displacement  thickness,  ft 

portion  of  drag  developed  up  to  any  particular  z  position,  lbs. 
fluid  kinematic  viscosity,  ft  "/sec 
fluid  density,  s!ugs/t>  * 
wall  skin-friction  coefficient 
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VIII.  5.  r.ibic  of  Nunifial  Solutions  of  Rise  Distance  liquation 
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